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Abstract 



In this work we develop a theory of Vessels. This object arises [LTJ [V] IBV] in the study of 
overdetermined 2D systems invariant in one of the variables, which are usually called time invariant. 
To each overdetermined time invariant 2D systems there is associated a vessel, which is a collection 
of system operators satisfying certain relations and vise versa. Such an invariance forces the theory 
of vessels to resemble a constant (classical) ID case NF B BC BGR, BFKD and as a result many 
notions are naturally redefined and most theorems are reproved in this setting. The notion of transfer 
function and its connection to the overdetermined 2D time invariant system (and the corresponding 
vessel) is one of the topics of this work. It is well known BC , BFKD that multiplicative structure of 
a transfer function of a ID system is closely connected to the decomposition of the state space into 
invariant subspaces of the state operator and we generalize this result to a wider class of functions. 
This class (denoted by I) arises as a class of transfer functions, which intertwine solutions of ODEs 
with spectral parameters CoLe:. At the end we present solution of factorization problems for finite 
dimensional case. 
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1 Introduction 



The theory of two-dimensional (2D) overdetermined time-invariant systems has been extensively devel- 
oped over the last 20 years; it is closely connected to the theory of commuting operators [LKMVj . |BV| . 
[V] , An overdetermined 2D continuous time- invariant linear i/s/o system is of the form 

( ^x(t u t 2 ) = A 1 x(t 1 ,t 2 ) + B 1 u(t 1 ,t 2 ) 
E: I £- x (t u t 2 ) = A 2 x(t u t 2 ) + B 2 u(t 1 ,t 2 ) 
{ y(ti,t 2 ) = Du(t u t 2 ) + Cx(t 1 ,t 2 ) 

where (u(ti, t 2 ), x(t\, t 2 ), y(t\, t 2 )) is the (input, state, output) triple and all the other symbols denote 
bounded operator on suitable Hubert spaces. Assuming continuously differentiable inputs, one obtains 
BV that the state space must be twice differentiable and enjoy the equality of mixed variables, from 
where algebraic relations are imposed on the operators of the system, and compatibility conditions 
(hence overdetermined system) on the input u(t\,t 2 ) and the output y(ti,t 2 ). More precisely, u(t\,t 2 ) 
and y(ti, t 2 ) satisfy algebraic equations (live on a curve) and an example of relation on the operators is 
commutativity of A\, A 2 . 

Using frequency domain analysis a notion of transfer function S^A,^) arises 

5(A) = D + C(XI-A 1 )- 1 B 1 

and the main question is how properties of such transfer functions and properties of the system operators 
(more precisely invariant subspaces of A\, A 2 ) are connected. Many standard structural properties, e.g, 
controllability, observability, minimality, pairing and adjoint system, cascade connection, equivalence 
and standard problems, e.g., pole placement, linear-quadratic-regulator problem H°° control for ID 
linear systems carry over for this setting. 

The study of time varying ID systems has produced a rich theory and the core of all these theories 
is the Sz.-Nagy-Foias theory of contractions on a Hilbert space [NFj . The analogue of transfer function 
in this setting is a lower block-triangular bounded operator. In [ABP], for example there is presented 
the development of a unified approach to time-varying dissipative linear systems, non-stationary Lax- 
Phillips scattering theory, and operator model theory for the infinite family of contractions. The abstract 
interpolation problem for the time-varying case is presented both in the de Branges-Rovnyak model 
formulation [dBRj, and in its coordinate-free, scattering-theoretic form, and its application to the time- 
varying version of the matrix right tangential Nevanlinna-Pick problem is studied. 

Probably, the simplest generalization of time varying ID systems is the study of 2D systems 
invariant in one direction. There are some works in this direction JLIJ IGa| in different settings. Our 
main inspiration comes from the article of M. Livsic [Li] and we actually continue this work. 

So, we suppose that our 2D systems are invariant in one of the variables (ti). We will use 
integrated form for of a system |La| and as a result there arises a continuous family of Hilbert spaces 
and semi-group acting between them. The invariance of the 2D system in one of the variables allows 
us to perform a partial separation of variables and to define a transfer function, depending on the 
corresponding spectral parameter (say A), which will, additionally, depend on the second variable (£2): 

S(X,t 2 ) = D(t 2 ) + C(t 2 )(XI - A 1 (t 2 ))~ 1 B 1 (t 2 ) 

A fundamental feature in the study of transfer functions and their factorizations is that algebraic equa- 
tions are now replaced by ODEs with a spectral parameter A. 
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The theory of these systems is interesting by itself, especially since it allows us to use frequency 
domain analysis in a time varying framework. It also has important connections with completely inte- 
grable nonlinear PDEs [Nj: the so-called Lax equation 

-7rAx{h) = A 1 (t 2 )A 2 (t 2 ) - A 2 (t 2 )A 1 (t 2 ) 
at 2 

appears naturally, and the passage from the input to the output ODE with a spectral parameter is 
analogous to the Backhand transformation. 

Let us introduce the detailed description of the contents of this work, highlighting the main results. 
Important notice is that we start the study of ii-invariant 2D systems from the more general integrated 
form, i.e., differential equations are presented as integral equations and as a result we obtain much 
weaker assumptions on the operators. Inspiration for doing it in this way comes from |Laj . As always 
the case in this passage, we consider evolution semi groups acting between continuous set of Hilbert 
spaces. This semi-group, after appropriate transformation of the Hilbert spaces and differentiation gives 
rise to operator A 2 (t 2 ) in section [3l 

There are seven sections in this work. After introduction, at the second section we introduce 
ovedetermined 2D systems and show how basic notions of system theory carry over. Following the 
ideas in [V] we introduce two notions (local and global) of approximate controllability and similarly two 
notions for observability and study their relations. On this basis, we study in section [3] guage quasi- 
similarity of (minimal) systems (vessels) . One of the interesting results of our work is that there exists 
a notion of differential vessel and we show that it is always possible to pass from integral to differential 
form and vice versa. Further, we present the notion of transfer function and its main properties. As a 
result we define a class X of intertwining functions, which is extensively studied at the last section El 

Next in section 0] we show that equivalence of transfer function for two vessels is equivalent to 
quasi-similarity between them. After that in section [5] the notion of adjoint system is presented. At 
section |6] we present basic operations on vessels: cascade connection, projection, compression, cascade 
decomposition. For the completeness of presentation we also discuss in section [7] Kalman decomposition 
in our setting that is very similar to the classical ID case. 



2 Over determined t\ invariant 2D systems 
2.1 2D systems invariant in one direction 

An overdetermined ^-invariant 2D system is a linear input-state-output (i/s/o) system, which consists 
of operators depending only on the variable t 2 \ in the most general case such a system is of the form 
P 

x(t 1 ,t 2 ) = e A ^*=)^- t ?)a ; (t?,t 2 ) + } 1 e^( t =)^-«)B 1 (t 2 ) U (y,t 2 )^ 

IYf • < * 2 

■ I x(h,t 2 ) ^ F(t 2 ,t°)x{h,t° 2 ) + J F(t 2 ,s)B 2 (s)u(t 1 ,s)ds 

> y{hM) = G{h)x{tx,h) + D(t 2 )u(ti,h) 
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where for Hilbert spaces £,£*,Ht 2 there are defined 

u(ti, t%) 6 £ - input, 
y{t\,t 2 ) e £* - output, 
x(t\,t 2 ) € Ht 2 - state, 

such that u(ti, t 2 ), y(ti, t 2 ) arc absolutely continuous functions of each variable when the other variable 
is fixed. The transition of the system will usually be considered from (t5,t§) to (tijtg)- Note that TCt 2 
are a priory different for each t 2 , and as a result F(t 2 ,t 2 ) has to be an evolution semi-group, i.e., it 
satisfies the following definition 

Definition 2.1 Given a collection of Hilbert spaces {Tit \t £ 1} for an interval ICR and a collection 
of bounded invertible operators F(s,t) : Ti s — » Ti t for each s,t £ I, we will say that F(s,t) is evolution 
semi-group if the following relations hold for all r, t, s G /: 

F(r,s)F(s,t)=F(r,t), 
F(t,t) = Id\ Ht . 

In order to be sure that all the formulas are meaningful we shall make the following regularity assump- 
tions 

Assumption 2.2 Internal regularity: 

1. Axfa) : Tit 2 — > Tit 2 , Bi (£2), B 2 (t 2 ) : £ — * T~Lt 2 , Cfo) ■ f~it 2 ~ y £* are bounded operators (for all t 2 ) 
and F(t 2 ,t 2 ) ■ Ti t o — ■+ Ti t2 * s & n evolution semi-group (see definition \2. 1}) . 

2. F(t 2 , s)B 2 (s) an d C(s)F(s,t 2 ) are absolutely continuous as functions of s (for almost all t 2 ) in 
the norm operator topology on £(£,Tit 2 ) and on C{£^,Ti t2 ) , respectively. 

Assumption 2.3 Feed through regularity: the operator D(t 2 ) : £ — > £* is an absolutely continuous 
function of t 2 . 

Since our entries u{t\,t 2 ) will be locally integrable as functions of t\, and as a result the first equation 
may be equivalently considered in the differential form, we shall work with systems of the following 
form: 

{^x(h,t 2 ) = A 1 (t 2 )x(t 1 ,t 2 ) + B 1 (t 2 )u(t 1 ,t 2 ) 
t 2 
x(t 1 ,t 2 )=F{t 2 ,4)x{t 1 ,4) + J F(t 2 ,s)B 2 (s)u(t u s)ds (2.1) 
y(h,t 2 ) = C(t 2 )x{h,t 2 ) + D(i 2 )u{h,t 2 ) 

2.2 Overdeterminedness and compatibility 

To ensure that the overdetermined systems equations (|2.I[) are compatible, we shall demand the equality 
of the two transitions for our system: 

f. (*?.*§) — » — » (*i,*2), 

2. (t°,t°) —>(ti,*!j) — >(ti,ta). 
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for arbitrary (t®,^), (ti,t 2 ). In the first case 



+ [ e A ^^-^B x {t 2 )u{p,t 2 )dp 



and in the second: 



ti 

+ F(t 2 ,4) [ e A ^ t ^B 1 (t° 2 )u(p 7 t° 2 )dp+ 



ta 

+ J F(t 2 ,s)B 2 {s)u(h,s)ds. 

The compatibility condition for free evolution (u = 0) results in 
or 

Alfa) = Ffa^A^Fitlh) (2.2.Lax) 

which is called the Lax equation [NJ and plays an important role in the theory of completely integrable 
non-linear PDEs. Note that it follows that the spectrum of Ai(t 2 ) is independent of t 2 . 

Inserting the Lax condition into x(ti,t 2 ) — x(tx,t 2 ) and rearranging the summands we obtain: 

ta ta 

e Ai(t 2 )( tl -t1) f ^( t2j s )fl 2 ( a ) u ( t j s ) ds _ /* F ( t2) a )fl 2 ( a ) u ( tl) s )d s = 



I a 

tl tl 



= F(i 2 ,t°) / e A ^ t ^B 1 (t° 2 )u(p,t Q 2 )dp~ / e^^-^Bxit^uip^dp 
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Multiplying this equality on the left by e W e reach 



*2 *2 

e~ M ^ J F(t 2 ,a)B 2 {a)u(t° 1 ,a)da - e~ M{t ^ J F(t 2) s)B 2 (s)u(h, s)ds = 

*2 *2 

tl tl 

= e -A 1 (t 2 )tx i? ( i2)f 0) y* e M{ti){ti- P ) Bi ^ U ( P ^ dp _ j e -Ai(t2)p Bl ( t2 } U ( P)t2 } dp 

Since u(p, a), e~ Al ^ t2 ^ p are absolutely continuous functions in each variable, we can rewrite it as an 
equality of iterated integrals of the derivatives (which are locally absolutely integrable as functions of 
one variable): 

/ / i [e- Al{t2)p F(t 2 , S )B 2 ( S )u(p, a)] dpda = JJf s [e- A ^F(t 2 , s)B 1 (s)u(p, a)] dsdp. 

*2 *1 ^1 ^2 

Notice that this integral equality is correct for all t1,t2,ti,t 2 . Moreover, the functions 
e~ Al ( t2 ^ p F(t 2 ,s)B 2 (s)u(p,s) and e~ Al ( t2 ' > ' p F(t 2 , s)Bi(s)u(p, a) are absolutely continuous and as a re- 
sult their derivatives are integrable functions, meaning that these two integrals are actually equal by 
Fubini's theorem to a two-dimensional integral over the rectangle {(p, a) t\ < p < t±, t\ < a < t 2 }. 
Thus it is equivalent to 

jL[e-MWr F fa a )Bi(8)u(p,s)] = ^[e- A ^F(t 2 , s)B 1 (s)u(p, a)} 

for almost all (p, a). This in turn is simply the following equation 

-A 1 (t 2 )e- A ^ p F(t 2 ,s)B 2 (a)u(p, a) + e ' A ^ p F{t 2l s)B 2 {a)^u{ Pl a) = 

op 

e~ A ^ p ^-[F{t 2 , a)B 1 (a)]u(p, a) + e~ A ^ p F(t 2l s )B 1 (s)^-u(p, a). 
Using (|2.2.Lax[ ) again we shall obtain 

-F(t 2 , 8)A 1 {8)B 2 {a)u(p, a) + F{t 2 , a)B 2 {s)-?-u(p, a) = 

op 

d d 
= -^[F(t 2 ,s)B 1 (s)]u(p, a) + F{t 2 , s)B 1 {s)—u(p, a), 

or, after multiplying on the left by F(a, t 2 ) and then substituting back the variables (p, s, t 2 ) — > (ii, ^2,^2) 

= 0. 
(2.3) 



B 2 {t 2 )—u(h,t 2 ) - B^—uihM) ~ (A!(t 2 )B 2 (t 2 ) + F(t 2 ,t° 2 ) — iFitlt^B^t^ufah) = 0. 
Oti ot 2 ot 2 



At this stage it is convenient to assume that we have factorization 
B 2 (t 2 ) = B(t 2 )a 2 (t 2 ), B 1 (t 2 ) = B{t 2 )a 1 (t 2 ), _ 

Ai(t 2 ) J B 2 (t 2 )+ J F(t 2j ^)^[i?(tS,t 2 ) J Bi(t 2 )]=-B(t a )7(i a ) 1 ' } 
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for some operators 

B{t 2 ) :£ -^H t2 , o- 2 (t 2 ), ai(t2),7(*2) ■£ ->f, 

where £ is another auxiliary Hilbert space. It is also important to postulate the following assumption 
in order to give meaning to the corresponding formulas: 

Assumption 2.4 External input regularity: 

1- 7(^2); ^2(^2) £ L[ oc (C(£ , £)) in the norm operator topology. 

2. <7i(t 2 ) S L(£,£) is absolutely continuous and invertible, in the norm operator topology. 

Expressed directly in terms of the operators B(t 2 ), cr i(t 2 ), o- 2 {t 2 ) with Bi(t 2 ), B 2 (t 2 ) eliminated, 
(12. 4[) becomes 

J-( J F(^,t 2 ) J B(t 2 )ai(t 2 ))+F(t°,t 2 )4 1 (t 2 )S(t 2 )<T 2 (t 2 ) + F (t°,t 2 )B(t 2 ) 7 (t 2 ) = 0. (2.5.0verD) 
at 2 

Then the condition (|2.3p becomes 

S(ta)[o-a(ta)- 5r «(ti 1 ia)-«ri(t a )ar«(ii,ta)+7(*2)]u(*i,i2) = 0. 

OTi OT 2 

A sufficient condition for this to hold (which is necessary in case B(t 2 ) is injective) is the input compat- 
ibility condition 

d d 

^{h)—u{t x ,h)-^i{h)^r «(*i s ta)+7(*aK*i.*a) = (2.6) 
oti ot 2 

Note that since u(t\,t 2 ) is the solution of PDE in the extended sense, it will be absolutely continuous 
as a function of t\ for almost all t 2l and conversely, it will be absolutely continuous as a function of t 2 
for almost all t\. 

The output y(ti, t 2 ) should satisfy the output compatibility condition of the same type as for the 
input compatibility condition (|2.6p . namely: 

d d 
(T2*(t2)g^y(ti,t 2 ) - au(t 2 )—y{t u t 2 ) + j*(t 2 )y(ti,t 2 ) = 0. (2.7) 

where similarly we have the following assumptions using an auxiliary output space £ * 

Assumption 2.5 External output regularity: 

1. 7*(i 2 ), o~2*{t 2 ) S L} oc (C(£*, £»)) m t/ie norm operator topology. 

2. <Ji*(t 2 ) £ £(£*,£*) is absolutely continuous and invertible, in the norm operator topology. 

So, inserting here t/(ii,t 2 ) = -D(t 2 )u(ii,t 2 ) + C(t 2 )x(t 2) t 2 ) we obtain that 

= K (£2)^-^*2)^+ 7* (*2)M*i,i 2 ) = 

= K(*2)^ - ^ + 7*(*2)][X>(t2)«(tl,«2) + C(t 2 )a:(t2,t2)] 

= ^(faJCte^sfei.ta) - ^i*(* 2 ) ^[CfeMia^)] +7*^2)^2)^2^2)+ 
+a 2 ,(i 2 )£>(t 2 )J r u(t 2 ,f 2 )- ( 7 1 4t 2 )J ? [ J D(t 2 )u(i 2l t 2 )]+7*(t 2 ) J D(t 2 )u(i 2l t 2 ). 
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Substituting here the first system equation from (|2.ip , we obtain (after omitting the notation of depen- 
dence on the variables) 

= <r 2 *C(A 1 x + Baxu) - oi*-^[Cx] + j*Cx + <r 2 *D-^u - a u [D'u + D-J^u] +j„Du = 
= ui*CA\x — ci, J^[Ci] -I- 7»Cx+ 

d 2 

' dti V1 *"dta 

Multiplying the second equation of the system (12. ip by C(t 2 ) and differentiating with respect to t 2 , one 

a 

hi , 



+ [(T2*D-M — a lm D-M- + (?2*CBoi - ai*D' + 7*D]u. 



can easily obtain that -M- [Cx] satisfies: 



ot 2 dt 2 
Inserting this into the last equation we obtain 

= <7 2 *CA lX - a u [CF]F- l x + CBo 2 u\ + ^Cx+ 

- auD-fe + o 2 *CBo x - a lx D> + 7 *D]u = 
= [a 2 *CA l - a 1 ^[CF]F- 1 + 7*C]x+^ { ' ' 

+[(?2*D - a^D-J^ + a 2 *CBa 1 - a^CBa 2 + a lt ,D' + j*D]u. 

The validity of this equation for the special case u(ti,t 2 ) = and an arbitrary initial x(ti,t 2 ) forces us 
to impose 

= o 2 *CA x F - au^-[CF] + j*CF. (2.9.0verD) 
ot 2 

With (|2.9.0vcrDj) in force, dUS} collapses to 

d d 

= \cr 2 *D- (Ji*D- h a 2 *CBoi - o\*CBo 2 + a^D' + ~/*D}u. 

oh ot 2 



On the other hand, u satisfies the input compatibility condition (|2.6| . so it is natural to assume that 
there is an operator D : £ * — > £* satisfying 

Assumption 2.6 External feed through regularity: the operator D(t 2 ) : £* — * is an absolutely 
continuous junction oft 2 , 

and which satisfies the following intertwining conditions, which will be called from now on the linkage 
conditions _ _ 

auD = Da l2 a 2 ,D = Da 2 , (2.10.Link) 

£>7 = cr 2 *CScri - a^CBa 2 + a^D' + 7 „L>. 

3 Vessel 
3.1 Definition 

Let us combine together all the formulas we have just developed. We define an (integral) vessel to be 
a collection of operator and spaces 

3W={A 1 {t 2 ),F{t 2 ,tl),B{t 2 ),C{t 2 ),D(t 2 ).D(t 2 )-a^ 
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satisfying regularity assumptions 12. 2[ 12. 31 12.41 12. 51 12.61 and the following vessel conditions: 

F(t 2 , tg)i4i(*jj) = A 1 (t 2 )F(t 2 ,t° 2 ) (12.2.Laxl) 
^(F(t°,t 2 )B(i 2 )<7 1 (t 2 )) + F(^,t 2 )A 1 (t 2 )B(i 2 ) < 7 2 (i 2 ) + F(t°,t 2 )B(t 2 ) 7 (t 2 ) -0 (l2.5.0verDI) 
^2*fe)Cfe)Ai(t 2 )F(t 2 ,^) -^(t 2 )^[C(t 2 )F(t 2 ,t°)] +7*(*2)C(t 2 )F(t 2 ,tE) = (12.9.0verDI) 

*j*D = Da h *2*D = Da 2 , fe.lO.Link t 
Z>7 = 02*CBo\ — o\*CBo 2 + (Ji„£) + 7*1?. 

It is naturally associated to the system IT (see (|2.1jl ) 

^x(tt,t 2 ) = Ai(t 2 ) x(t u t 2 ) + B(t 2 )ai(t 2 ) u(h,t 2 ) 

I^- ( x{h,t 2 ) = F(t 2 ,4)x(t u 4) + f F(t 2 ,s)B(s)a 2 (s)u(t 1>S )ds (3.1) 

t° 

, y(h, h) = C{t 2 ) x(h, t 2 ) + D{t 2 )u{ti,t 2 ). 

with absolutely continuous inputs and outputs, satisfying compatibility conditions (|2.6|) . (|2.7|) for almost 
all (ti,*a): 

o-2(t2)TOr«(ti,t2) - CTi(* 2 )^u(ti, t a ) + 7(*2)«(*i, *a) = 0, 
ct 2 * (t^^yih , t 2 ) - a u (t 2 ) (tj , t 2 ) + 7* (t 2 )tf(*i , *a) = 0. 

We shall further name these conditions as follows. (|2.2.Lax[ ) is the Lax equation. (|2.5.0vcrD[) and 
(|2.9.0verD|) are the input and the output vessel conditions, respectively. (|2.10.Link[) is the linkage 
condition. 

3.2 Gauge quasi-similarity of vessels 

As in the classical case, in order to deal with some classification of systems (to be defined later) we 
need a notion of minimal systems. In the ID case there is only one natural notion of approximate 
controllability and observability. In the case of 2D t\ invariant systems, on the other hand, there are at 
least the following two notions 

Definition 3.1 System IT Ii3.1\) is called locally approximately controllable at t 2 if 

C t2 = {h e H t2 I 3*i e R, (u, x,y)&T : x(0, t 2 ) = 0, x(h,t 2 ) = h} 

is dense in 7it 2 ■ 

Here T stands for the set of system trajectories (u, x, y) of (|3.1[) with compatibility ODEs l^TT)) and 
with the spectral parameter A. 

Definition 3.2 System IT \3.1\) is called approximately controllable at t 2 if 

C t2 ={he H t2 I 3*i e R,{u,x,y) e T : a;(0, 0) = 0, x(t u t 2 ) = h] 

is dense in TL t2 ■ 
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Similarly to the classical case one immediately obtains that 

C i2 =span{ImAi(f 2 )B(t 2 ) | j = 0, 1, 2, . . .}, 

C t2 =spm{ImF(t 2 ,s)A{(s)B(s) \ j = 0, 1,2, . . . , s e R}. 

Notice that -F(i2, t' 2 )Ct> — Ct 2 and since -F(i2j ^) are bounded invertible operators, the density of Ct 2 for 
any value of t 2 implies density for all t 2 , which means that the notion of approximate controllability is 
independent of t 2 . 

Since F(t 2 ,t 2 ) — Id, we obtain that Ct 2 Q Ct 2 for all t 2 and consequently, if Ct 2 is dense in Ht 2 
then so is Ct 2 . But actually the converse also holds 

Theorem 3.1 For the system IS hS.l\) , satisfying regularity assumvtions \2.2[ \2.!A\2.4\ \2.5[ \2.b\ the local 
approximate controllability for arbitrary t 2 is equivalent to approximate controllability for all t 2 ■ 

Proof: We have seen that 

= 1~L t o => Cfo = H t o C t2 = H t2 

and we need to prove the converse. Suppose that for a fixed t 2 we have approximate controllability, 
which means that 

V F(t° 2 ,s)A{(s)B(s)e = H t » 

n>Q,e££ 

or what is equivalent 

B*(s)F*(4, s)Al(t° 2 )h = 0, Vfte H t o. 
In the same manner local approximate controllability means that 

V A?(4)B(t° 2 )e = H t o B*(t°)A™(4)h = 0, Vfte H t o. 

n>0,ee£ 

Consider now a function with values in £ for each h 6 H+o 

2 

u h {\t 2 ,4) = B*(t 2 )F*(t° 2 ,t 2 )(XI - Alitl))- 1 ^ u h (X,t 2 ,t° 2 ) = B*(t° 2 )(\I - AUtl))-^ 

which is analytic at the neighborhood of A = oo. This function satisfies the following differential equation 
(adjoint of the input vessel condition (|2.5.QverDp ) 

&u h (\ t t 2 ,4) = -^-[B* (t 2 )F* (4,t 2 )(\I - A\(t 2 \))~ 1 h] = 

= a^ 1 [-a 2 B*(t 2 )F*(4,t 2 )A* 1 n (4)(XI - At(4))- 1 h + 7 B*(h)F*(4,t 2 )At n (4)(XI - A^))-^} = 
= a^[-a 2 B*(t 2 )F*(4, t 2 )(At n (4) - XI + XI)(XI - A\{4))- l h + ~fu h (X, t 2 ,t° 2 )} = 
= a^l-asX + j]u h (X, t 2 ,4) + a^a 2 B* (t 2 )F*(4, t 2 )h. 

denoting by §(X,t 2 ,4) the fundamental solution of the differential equation 

-^U(A, t 2 , t° 2 ) = a? [-* 2 X + 7 ]$(A, t 2 , 4) 
at 2 

and using variation of parameters, we shall obtain that 

u h (X, t 2 ,4) = $(A, t 2 ,4)[$ $-!(A, y, 4)a^a 2 B*(y)F*(4, y)h + u h (X, 4,4)] = 
= ®{X,t 2 ,4)[f t ^-\X,y,4)a- 1 a 2 B*\y)F*^^^ 

= $(A, t 2 ,4) f| <J> _1 (A, y, 4)a^a 2 B* (y)F* (4, y)h + $(A, t 2 ,4)B* (t 2 )(XI - A 1 }(4))~ 1 h, 
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which is a sum of an analytic (in A) £- valued function and of a function with poles (for A in the spectrum 
ofAl(t° 2 )) 

$(\,t 2 ,4)B*(t° 2 )(\I - AUtl))-^. 

So, if there is a vector h G H t o such that Uh(X,t2,t 2 ) — f° r each value of A, it has to vanish the part 
containing the pole first. This means that for all A out of the spectrum of A\(t 2 ) 

<S>(X,t 2 ,t° 2 )B*(t° 2 )(XI-A* 1 (t° 2 ))- 1 h = 

and since &(X,t 2 ,t 2 ) i s invertible we obtain, 

B*(t° 2 )Al n (t° 2 )h = 0, 

which is local approximate controllability. □ 
The following notions are natural generalizations of the ID case: 

Definition 3.3 The system IT is locally observable at t 2 if 

0£ = {h e H t2 | x(0, t a ) = h, («, x, y) G T, u{h,t 2 ) = => y{h,t 2 ) = for all t x ) = {0}, 

and the system is observable at t 2 if 

&i 2 ={heU \ x(0,0) = h, (u,x,y) G T, u=0^y=0 for all h} = {0}. 

Again, similarly to the classical case one obtains that 

Ok =n„ eN Ker(7(t 2 )Ay(t 2 ), 

Oh = n„ eM Ker<7( S )^( S )F(M 2 ). 

From the property of the evolution semi-group F(t 2 ,t 2 ) = I we obtain that 0^ 2 2 and that 
local observability at t 2 (i.e., 0^ = {0}) implies observability at t 2 (i.e., 0^ = {0]}). Additionally, 
Oh. = Ob F(t' 2 ,t 2 ) implies that observability is independent of t 2 , i.e., if for any t 2 , Om — {0} then 
0^ = {0} for all t 2 . Finally, a parallel to the controllability theorem. 

Theorem 3.2 For the system IT iS.l)) , satisfying regularity assumvtions \2.2[ \2.!A\2.4\ \2.5l \2.b\ the local 
observability for arbitrary t 2 is equivalent to observability for all t 2 . 

Proof: We have seen that 

ojo = {0} =► 6\ = {0} => 6i 2 = {0} 

and it is remained to show the converse. Similarly to the proof of local approximate controllability for 
each t 2 , supposing observability at a fixed t 2 

Oi 2 = P| KevC(s)A^(s)F(s,t 2 ) = {0}, 

nGN.sGR 

we may consider a function of A with values in £ 

y h (X, t 2 ,t° 2 ) = C(t 2 )F(t 2 ,t a 2 )(XI - A^))- 1 ^ 
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and the output vessel condition (|2.9.0vcrD[) 

V2*(t2)y Al (t°)h(t2) - vi*(t 2 )-^y h {t 2 ) + 7*(t 2 )i/ h (t 2 ) = 0, y h (X, t° 2 , t° 2 ) = C(t° 2 )(\I - A 1 {t° 2 ))- 1 h 

and follow the same lines as in the proof of local approximate controllability. □ 
The system JS is called minimal if it is both approximately controllable and observable, i.e., 
(C t2 = O t2 =H). 

A natural notion of similarity arises, which is used to classify (usually minimal) systems and 
corresponding vessels. Two vessels 

32J = {A x (t 2 ), F(t 2 , 4), B{t 2 ), C(t 2 ), D(t 2 ), D{t 2 ); a x {t 2 ), a 2 (t 2 ), 7 (i 2 ), tr u (t 2 ), ct 2 .(* 2 )7. Wfa , £, £*, £ £.) 
323 = (Ii(t 2 ), F(t a , * 2 ), B(t 2 ), C(t 2 ), D(t 2 ), b{t 2 );a x {t 2 ), a 2 (t 2 ), 7 (i 2 ), a u (jb 3 ), a 2 *(t 2 )j*(t 2 );iit 2 ,£, £*,£, £*) 

are called gauge quasi-similar (or kinematically quasi- similar), if there exists a (possibly unbounded) 
linear operator T(t 2 ) : D(T(t 2 )) — > Wt 2 , with a dense domain D(T(t 2 )) C 7i t2 , which is 1-1, with dense 
range, and satisfies the following intertwining conditions 

iii(t 2 )T(i 2 ) =T(t 2 )Ai(t 2 ), 
F(t 2 ,tO)T(t°) = T(t 2 )F(i 2 ,i°), 

i(*2) =T(t 2 )B(t 2 ), 
C(t 2 )T(t 2 ) = C{t 2 ). 

Moreover, we shall demand that 

4 C D(T(t 2 )), F(f a ,t°)£>(T(f°)) C D(T(*a)). (3.3) 

These conditions arc necessary in order to obtain reasonable definitions in (|3.2[) . For example, F(t 2 , i 2 )F(i°) = 
T(t 2 )F(t 2 ,4) requires F(h,t%)D(T(l%)) C D(T(t 2 )). The first and the third conditions require 
ImB(t 2 ),ImAi(i 2 ) C D(T t2 ), for which it is enough to demand C 4a C D(T(t 2 )). 

When it is the case that D(T{t 2 )) = Tit 2 for all t 2 , and T(t 2 ) is everywhere defined bounded and 
onto (with bounded inverse as a result), we say that the vessels 323, 323 are similar. 
Remark: 1. For the finite dimensional case dim£ < oo, as in the classical case the notions of similarity 
and of quasi-similarity coincide. 

2. Given a vessel 323 and a family of invertible bounded operators T(t 2 ) : Tit 2 — ► Wt 2 the formulas (13. 2|) 
define a new vessel 353 that is gauge similar to 353. 

3.3 Differential vessels 

Suppose that the system IT, defined by (|3.1[) has identical inner spaces Ht 2 for all t 2 , i.e., Ht 2 = H. 
Suppose also that the evolution semi-group F(t 2 , t§) is absolutely continuous (for a fixed t 2 ) in the norm 
operator topology of B(H,H), then its generator A 2 (t 2 ) is for almost all t 2 a bounded operator and 
satisfies 

A 2 (t 2 ) = ^-F(t 2 ,t° 2 )F- 1 (t 2 ,t° 2 ). 
at 2 

Rewriting all the vessel conditions in the differential form using j4 2 (t 2 ), we obtain a differential vessel 
2)23 = (Ax(t 2 ), A 2 (t 2 ), B(t 2 ), C(t 2 ), D(t 2 ), D(t 2 );ai(t 2 ),a 2 (t 2 ), j(t 2 ),a u {t 2 ), a 2 *(t 2 )^{t 2 );H, £, £*,£, £*) 
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which satisfy the following axioms: 



^Ai(i a ) = A 2 {t 2 )A x {t 2 ) - A 1 (t 2 )A 2 (t 2 ) 

^(S(t 2 ) CTl (f 2 )) - A 2 {t 2 )B{t 2 )a l (t 2 ) + i4i(t 2 )B(*2>2(t 2 ) + B{t2)l{t 2 ) = 
^(^(t 2 )C(* 2 )) +vi*{t2)C{t 2 )A 2 (t 2 )+a 2 *{t 2 )C{t 2 )A 1 [t 2 ) + 7 *(i 2 )C(i 2 ) = 



ui*Z) = Da\, a 2 *D — Da 2 , 



Dj = o^CBox - a^CBa 2 + a u D' + 7*L>. 



and the following regularity assumptions, which are obtained from the assumptions 12. 2[ 12.31 12.41 12. 51 
12.61 by differentiating" 

Assumption 3.4 1. Internal regularity: Ai(t 2 ), A 2 {t 2 ) : Ttt 2 — > 7it 2 , Bi(t 2 ), B 2 (t 2 ) : £ —> 7it 2 , 
C(i 2 ) : Wt 2 — > £* are bounded operators (for all t 2 ), 

2. Feed through regularity: the operator D(t 2 ) : £ — ► £*, D(t 2 ) :£*—>£* are absolutely continuous 
function of t 2 , 

3. External regularity \2.4\2.5\ 

The differential vessel is associated with the system 



This motivates the following definition 

Definition 3.5 Vessel 323 will be called differentiable if Ht 2 = H for all t 2 and i* 1 ^,^) is an abso- 
lutely continuous function in the norm operator topology of B(Ti,,Ti,). 

Then the following proposition holds 

Proposition 3.3 Any vessel 327 is gauge similar to a differentiable vessel. 

Proof: Let us take T(t 2 ) = F(t%,t 2 ) : H t2 -> H t o, which is an absolutely continous, and let us find out 
what kind of a vessel is obtained. First notice that T(t 2 ) are all mappings onto the same space "H t o, 
which we denote by H. t o = TL. Further, we see that 



£-x(t u t 2 ) = Atit^xit^h) + B(t 2 ) <n(t 2 ) u(t u t 2 ) 
£>£: < ^Lx(h,t 2 ) = A 2 (t 2 )x(t!,t 2 ) + B(t 2 ) a 2 {t 2 ) u(h,t 2 ) 
^ y{ti,t 2 ) = D(t 2 )u(h,t 2 ) + C(t 2 )x(h,t 2 ) 



(3.4) 



and compatibility conditions for the input/output signals: 



o-2(t 2 )^-u(t!,t 2 ) - ax(t 2 ) -^u(h,t 2 ) + 7(t 2 )u(ii,t 2 ) = 0, 
<J2* (h)^-y(h , t 3 ) - o-u (t 2 ) ^y (h , t 2 ) + 7 * (t 2 )y(h ,t 2 )=0. 
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Thus we obtain a differcntiablc vessel 



S2J = (F(t° 2 , t 2 )A 1 (t 2 )F(t 2 , t° 2 ), I, F(t° 2 , t 2 )B(t 2 ), C(t 2 )F(t 2 ,t 2 ),D(t 2 ),D(t 2 y y 
°i(t 2 ), &2(t 2 ), j(t 2 ), v^fo), cr 2 *(t 2 )^(t 2 );H, £,£*,£,£*) 

with trivial evolution semi-group /, which is obviously absolutely continuous. □ 

Notice that this vessel is of a very special form. Its evolution semi-group is trivial and as a result, 
differentiating this vessel, we shall obtain that the generator of the semi group A 2 {t 2 ) = is trivial. 
Consequently, from the Lax equation, Ai(t 2 ) = A\ becomes a constant operator. 

Let us for the completeness of presentation explicitly write down the notion of quasi similarity 
of two differential vessels. Gcvcn two vessels 

J)W=(Ai{t2),A i (h),B{t 2 ),C{fr),D^ 

and 

2)9J= (Ii(t 2 ),i 2 (t 2 ),i?(t 2 ),<5(i 2 ),i}^ 

we will say that they are quasi-similar if there exists a (possibly unbounded) linear operator T(t 2 ) : 
D(T(t 2 )) — > 7it 2 , with a dense domain D(T(t 2 )) C Ti. t2 , which is 1-1, with dense range, absolutely 
continuous (in the norm operator topology) and satisfies the following intertwining conditions 



f ii(t 2 )T(t 2 ) 
A 2 {t 2 )T{t 2 ) 

I C(t 2 )T(t 2 ) 



= T(t 2 )A 1 (t 2 ) 
= T(t 2 )A 2 (t 2 ) 

= T(t 2 )B(t 2 ), 
= C{t 2 ). 



(3.5) 



We shall also demand that 



C t2 C D(T(t 2 )), 



lmA 2 (t 2 ) C D{T{t 2 )). 



If T(t 2 ) is an invertible bounded operator, we shall say that two differential vessels are similar. 



(3.6) 



3.4 Separation of variables and the notion of transfer function 

One of the reasons why overdetermined systems invariant in one direction are interesting is the possiblity 
tp perform a partial separation of variables. Taking all the trajectory data in the form 

u(tut 2 ) = u x (t 2 )e xt \ 
x(t u t 2 ) =x x {t 2 )e xt \ 
y(h,t 2 ) =y x {t 2 )e xt \ 

we arrive at the notion of a transfer function. Note that u(t\,t 2 ), y{t\,t 2 ) satisfy PDEs, but u\(t 2 ) 7 y\(t 2 ) 
are solutions of ODEs with a spectral parameter A, 



d 

\a 2 (t 2 )u\(t 2 ) - <Ji(t 2 )—u x (t 2 ) + j(t 2 )u x (t 2 ) 



0. 



Xv2*(t 2 )y\(t 2 ) - (Ti*{t 2 )—y x {t 2 ) + l*{t 2 )y\{t 2 ) = 0. 



(3.7) 
(3.8) 
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The corresponding i/s/o system becomes 

x\fo) = (XI - A.ih^-^i^a^uxih) 
x(ti,t 2 ) = F(t 2 ,T 2 )x(ti,T 2 ) + J F(t 2 ,s)B(s)a 2 (s)u(ti,s)ds 
V\{t 2 ) = D(t 2 )u x (t 2 ) + C{t 2 )xl{t 2 ) 
The output y\(t 2 ) = D(t 2 )u\(t 2 ) + C(t 2 )x\(t 2 ) may be found from the first i/s/o equation: 

y\(t 2 ) = S(X,t 2 )u x {t 2 ), 

using the transfer function 

S(X,t 2 ) = D(t 2 ) + C(t 2 )(XI - Aitejr^teMia). (3.9) 

Here A is outside the spectrum of Ai(t 2 ), which is independent of t 2 by (|2.2.Lax[) . We emphasize that 
S^A, t 2 ) is a function of t 2 for each A (which is a frequency variable corresponding to t\). 

Theorem 3.4 A transfer function S(X,t 2 ), defined by i3.9\) has the following properties: 

1. For almost all t 2 , S(X,t 2 ) is an analytic function of X in the neighborhood of oo, where it satisfies: 

S(oo,t 2 ) = D(t 2 ) 



2. For all X, S(X,t 2 ) is an absolutely continuous function oft 2 . 

3. For each fixed X, multiplication by S(X,t 2 ) maps solutions of 
Xa 2 {t 2 )u — o~i(t 2 )j^ +j{t 2 )u — to solutions of 

X(j 2 (t 2 )y - o-i (i 2 )^ + l*{t 2 )y = 0. 
Proof: Those are easily checked properties, following from the definition of S(X,t 2 ): 

S(X,t 2 ) = D(t 2 ) + C(t 2 )(XI - Ai(t 2 ))~ 1 B(t 2 )o-i{t 2 ). 

When A — > oo, since all the operators are bounded the second summand vanishes and we obtain 
S(oo,t 2 ) = D(t 2 ). Moreover, it will be an analytic function of A, when A > ||Ai(t2)|| and we obtain the 
first property. 

In order to understand the second property let us rewrite S*(A, t 2 ), using the Lax equation in the 
following way: 

S(X,t 2 ) =D(t 2 )+C(t 2 ){XI-A 1 (t 2 ))- 1 B{t 2 )ai(t 2 ) = 

= D(t 2 ) + C(t 2 )(XI - F(t 2 ,tO)^ 1 (tO)F(tO,t2))- 1 -B(t 2 )(Ti(t2) = 
= D{t 2 ) + C{t 2 )F{t 2 ,t%){XI - A^t^Fitl t2)B(t 3 )<T 1 {t 2 ). 

The functions C(t 2 )F(t 2 , t 2 ), F(t 2 ,t 2 )B(t 2 ), <Ji(t 2 ) are absolutely continuous in appropriate spaces, thus 
their multiplication too and we obtain the second property. 

The third property is a direct result of our construction. □ 

Remark: for the case dim7i t2 < oo, we obtain that S(X,t 2 ) is a rational (off the spectrum of Ai(t 2 )) 
in A function for all t 2 . 
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3.5 Class X of intertwining functions 

We saw in previous section (thorem 13. 4p that transfer functions are very natural objects to study and 
have three important properties. Suppose that we start from two ODEs the input [3~7l and the output 
13.81 and denote the fundamental solution for them by Q{X,t 2 ,t^) an d by ^(A, t 2 , t 2 ) respectively, then 

S{\, i 2 )$(A, t2,4) = $ *(A, t 2 ,T 2 )S(\, t°) (3.10) 
and S(X,t 2 ) satisfies the following ODE 

J-S(A, h) = a^ 1 {t 2 ){a 2 , (t 2 )X + 7* (t 2 ))S{X, t 2 ) - S{\, t 2 )a^ {t 2 ){a 2 (t 2 )X + 7 (t 2 )). (3.11) 
ot 2 

We recall [CoLe] that from the fundamental theory of linear differential equations that for each equation 
there correspond and invertible matrix (or operator) function <f>(t 2l t 2 ) which obtains value / for a fixed 
value of t 2 — t 2 and any other solution u(t 2 ), satisfying u(t 2 ) = uq is just of the form 

u{t 2 ) = (j>{t 2 ,t° 2 )u . 

Some of the simple properties of the fundamental matrix <I>(A, £2,^2) are as follows. Notice that 
<I>(A, £2, ^2) can be replaced by $*(A, £2, ^2) with a corresponding change of operators. 



Lemma 3.5 The fundamental matrix ^(A, t 2l t 2 ) in \3. 9f ) satisfies: 

J-$(A, t 2 ,t° 2 ) = a^(Xa 2 (t 2 ) + 7 (t2))*(A, t 2 , t°), (3.12) 

^-$~\X,t 2 ,t° 2 ) = -<p- 1 (X,t 2 ,t a 2 )a^ 1 (Xa 2 (t 2 )+ 1 (t 2 )), (3.13) 

^*(X,t 2 ,t a 2 ) = <P^X,t 2 ,t° 2 )CXaUt 2 ) + l*(t2))^ 1 *, (3.14) 
ot 2 

^- 1 %XMA) = HM{t2) + i*{t 2 ))^ 1 *Hx,t 2 Ar 1 \ (3-i5) 

ot 2 

Proof: Conjugating (|3.7p and using a formula for the derivative of the inverse. □ 



Since we will intensively work with such functions, we define their class X as follows 
Definition 3.6 The class 

I = Z((7i,0'2,7; cr l*, '2*,7*) 

is a class of functions S(X,t 2 ) of two variables, which are 

1. analytic in the neighborhood of X = 00 for all t 2 , 

2. absolutely continuous as functions of t 2 for almost all X, 

3. map solutions of the input ODE \3. 7| ) with spectral parameter X to the output ODE \3.8)) with the 
same spectral parameter (i.e., they satisfy the equation (3.10\) ) 
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4 Main theorem of gauge quasi similarity 

The following result is an analogue in our framework of the standard quasi similarity theorem for minimal 
systems piBC]. 

Theorem 4.1 Assume that we are given two minimal (integral) vessels 353, 353 

353 = (Ai(t 2 ), F(t 2 , t° 2 ), J3(t 2 ), C(t 2 ), D(t 2 ), 5(t 2 ); ffi(t 2 ), <T 2 (t 2 ), l(t2),o-u{t 2 ), <T 2 *(t 2 ) 7 *(i 2 );ft t2 , £, £*,£, £*) 
3W=(A 1 (t 2 ),F{t 2 ,t° 2 ),B{t 2 ),C{t 2 ),D{t 2 ),b(t 2 ^ 

with transfer functions S(X,t 2 ), S(X,t 2 ). Then the vessels are gauge quasi similar iff S(X,t 2 ) = S(X,t 2 ) 
in a neighborhood of X = oo. 

Proof: The easy direction of the statement considers the case when there exists T(t 2 ) : Ht 2 — * ~Ht 2 , 
responsible for quasi-similarity, i.e., satisfying (|3.2[) and (|3.3|) , then 

S(X,t 2 ) = D(t 2 ) + C(t 2 )(XI - A 1 (t 2 ))- 1 B(t 2 )a 1 (t 2 ) = D(t 2 ) + C(t 2 )T(t 2 )(XI - Aifa))^ 1 B(t 2 )(Ti(t 2 ) = 
= D(t 2 ) + C(t 2 )(XI - A 1 {t 2 ))- 1 T{t 2 )B{t 2 )a 1 {t 2 ) = D(t 2 ) + C(t 2 )(XI - A^))- 1 B^a^) = 
= £(A,t 2 ). 

Fo the converse, we obtain first that the values at infinity of the two functions are equal: 

D{t 2 ) = D{t 2 ) => D(t 2 ) = D(t 2 ), since o"i,cri* are invertible. 

By looking at the Taylor coefficients in the power series expansions of functions S(X,t 2 ), S(X,t 2 ) at 
infinity, we obtain 

C{t 2 )A^(t 2 )B{t 2 ) = C{t 2 )A r l(t 2 )B{t 2 ), for all n = 0, 1, 2, . . . (4.1) 
Since 23 is approximately controllable, the set 

N 

C t2 = {^F(t 2 ,s0^i 1 (^)B(s J )e I | ni ,N G N,e< G £,s, G K} 
i=i 

is dense in Ti. t2 . Define next T : Ct 2 — > Ti.t 2 by 

JV N „ 

i=l i=l 
Then T is obviously a linear transformation, provided it is well defined, i.e. we have to check that 

N N _ 

= F{t 2 , Sl )A n ^ => = J] F(t 2 , s;)!^ 

i=i i=i 

Since by assumption 53 is observable, to show = Yli=i F{t 2 , s^A^ {si)B{si)ei is the same as showing 
that 

C(t 2 )A^(t 2 )(J2 F{t 2 , Sl )A n ^ (8i)6(si)ei) = 0, for all k = 0, 1, 2, . . . 

j=i 

This is done with the help of the following lemma 
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Lemma 4.2 For each t 2 ,s the following equality holds: 

C(t 2 )(XI - A 1 (t 2 ))- 1 F(t 2 ,s)B(s)a 1 (s)e = C{t 2 ){XI - A 1 {t 2 ))- l F{t 2 ,s)^{s)o- 1 {s)e 
Proof: Simple calculations using the first vessel condition 12.5 . vcrDl result in 
§- 8 {C{t 2 ){XI - A 1 {t 2 ))- l F{t 2 , s)B{8)a 1 (s)) = 

= -C(t 2 )(XI ~ A 1 {t 2 ))- 1 F{t 2 , s)B{s)[a 2 {s)X + 7 (s)] + C(t 2 )F(t 2 , s)B(s)a 2 (s), 
Remember that the function $(A, t 2 , ^2) satisfies (|3.13p . so using variation of coefficients 
C(t 2 )(XI - Ai(t 2 ))~ 1 F(t 2 , s)B(s)ai(s) = K(X, t 2 , s)^ 1 (X, t 2 , s) 

where 

^- s K(X,t 2 ,s) = C(t 2 )F(t 2 ,s)B( S )a 2 (s)<S>(X,t 2 ,s). 

Particularly, K(X 7 t 2l s) is 

K(X, t 2 , s) = K(X, i°, s) + [' ^-K(X, t 2 , s)ds = 

= K(X, 4, s)+ f C(t 2 )F(t 2) s)B(s)o- 2 (s)$(X, t 2) s)ds 

and all its poles are at the first term K(X,t^,s). The same considerations, applied to C(t 2 )(XI — 
A-i(t 2 ))~ 1 F(t 2 , s)B(s)ai(s) results in a function K(X, t 2 ,s), whose poles are also at the term K(X, t®,s): 

K(X, t 2) s) = K(X, t° 2) s)+ r C(t 2 )F(t 2 , s )B( s )o- 2 ( s )$(X, t 2 ,s)ds. 

Jt% 

Then 

C(t 2 )(XI - A 1 (t 2 ))~ 1 F(t 2 , s)B(s)ai(s) - C(t 2 )(XI - Ai(t 2 )) _1 F(<2, s)_B(s)cti(s) = 
- lK(X,t 2l s)-k(X,t 2 ,s)}<t>- 1 (X,t 2l s) = 

= [K(X,t 2l t 2 ) + J t l^K(X,t 2 ,y)dy - K(X,t 2l t 2 ) - j; 2 ^k(X,t 2 ,y)dy]<p-\X,t 2 ,s) = 
= [Jl JLK(X,t 2 ,y)dy - g ^k(X,t 2 ,y)dy]^-\X,t 2 ,s), 

in other words, this difference is an entire function of A. On the other hand, two functions 

C(t 2 )(XI - Axfo))^ 1 F(t 2 , s)B(s)cr 1 (s), C(t 2 )(XI - Aifo))^ 1 F(t 2 , s)B(s)ai(s) 

are zero at infinity (i.e., globally bounded). Taking a paring with an arbitrary linear functional and the 
operator applied to arbitrary vector, by Liouville's theorem this difference is constant and is equal to 
the value at infinity for each such pairing. Thus the operator itself is zero. □ 

One of the consequences of this theorem is that Taylor coefficients around infinity of the functions 
C(t 2 )(XI - Aifo))^ 1 F(t 2 , s)B(s)ai(s), C{t 2 )(XI - A 1 {t 2 ))~ x F(t 2 , s)&(s)a 1 {s) 
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are equal, and consequently, 

C{t 2 )A 1 {t 2 ) n F{t 2 ,s)B{s)a 1 {s) = C{t 2 )A 1 {t 2 ) n F{t 2 ,s)h{s)a 1 {s) 
Using this result, the fact that T{t 2 ) is well defined is immediate, because 

C(t 2 )AUt2)(EtiF(h, St )A n 1 '( Sl )]3( Sl )e l ) = £ C(t 2 )F(t 2 , Si )A? +k = 

i=0 

N „ N „ 

- E C(t 2 )^(t 2 , S ^i l+fc (^)£?(^)e* = C(t 2 )Af(i 2 )(X; F(t 2 ,s t )A^(s,)B( Sl )e t = C(t 2 )A k (t 2 )0 = 0. 

i=0 i=0 

One checks that T is one-to-one by using the observability of 03: 

KerT = n„ eN ,.€R Ker C( S )A?(s)F(s, t 2 ) = 6£, 

which is trivial by assumption of minimality. That T has dense range follows from the approximate 
controllability of 2J. This finishes the proof. □ 

Since the notions of similarity and quasi-similarity are identical for finite dimensional vessels, we 

obtain 

Corollary 4.3 Two finite dimensional vessels (with rational in A transfer functions) are gauge similar, 
if and only if, the transfer functions are identical at the neighborhood of A = oo. 

5 Adjoint system 

The notion of the adjoint system is very useful in system theory. It is obtained from a simple observation 
that applying adjoint to the vessel conditions gives rise to a new set of conditions on adjoint operators, 
which are almost vessel conditions. Moreover (£*)* is actually S by a trivial change of coordinates 
x — > —x on the state space. Here is the precise definition. Given a system in the differential form 

^x{t u t 2 ) = A^xih,^) + B(t 2 ) (Ji(t 2 ) u(h,t 2 ) 
3S:< ( 5£*(*i»*2) = A 2 (t 2 )x(h,t 2 ) + B(t 2 ) a 2 (t 2 ) u(h,t 2 ) (5.1) 
y(t u t 2 ) = D(t 2 )u(ti,t2) + C(t 2 )a;(ti,t2) 

and associated vessel 

W={A 1 {t 2 ),A 2 {t 2 ),B{t 2 ),C{t 2 ),D{t 2 ),D{t 2 )-,a 1 {t^ 
it is natural to introduce the adjoint system: 



DT,* : { 



-^x,{t 1 ,t 2 )=A* 1 (t 2 ) x*{t u t 2 ) + C*{t 2 )at*u*(t u t 2 ) 

-^x,{t u t 2 ) = A* 2 {t 2 ) x*{t u t 2 ) + C*{t 2 )a^u,(t u t 2 ) (5.2) 



. y*(t u t 2 ) = B* x*(t u t 2 ) + D*(t 2 ) u.(*i, t 2 ), 
which is associated to the vessel 53* given by 

2?* = (--A*) -A* 2 , -C* , B* , D* , D* ; cr**, <7 2jk , -7* - — cr*„, cr*, crj, -7* - ^^l^' £*,£*,£*,£), 
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where all the operators are functions of t 2 and satisfy the following axioms: 
JL A* — A* A* — A* A* 

-iS C yt*) - A* 2 C*a* u + A$C*a^ + C*( 7 : + ^trf.) = 
^(*tB*) - a* x B*A\ + a* 2 BA\ + 7 *B* = 
a* x D* =D*a* u , a* 2 D* = D*a* 2 , 

D *(-~/t - J^i*) - -P$B*C*a* u + oiB*C*o& - al^-D* - ( 7 * + &o?)D*. 
Moreover, the transfer function of the adjoint vessel S*(p,,t 2 ) maps solutions of the adjoint input ODE 

dt 2 dt 2 

with the spectral parameter \x to solutions of the adjoint output ODE 

[fT2> ~ a *~k 7 * ~ Jr 2 al]y * {t2) = (5 - 4) 

with the same spectral parameter. In the language of fundamental matrices, if one denotes the funda- 
mental solution for u*(t 2 ), y*{t 2 ) as &((j,,t 2 ,t 2 ) an d t 2 , t 2 ), respectively, then (similarly to (|3.10p ) 

S*(fJ,,t 2 )^(iJ,,t 2 ,T 2 ) = ^*(^,t 2 ,T 2 )S^(p,T 2 ) 

Notice also that the equation for S*(X,t 2 ) similar to (|3.9p is: 

S*((x, t 2 ) =D*- B*(fxl + AD^CTa^ 

And as in the case of constant operators, one obtains that 2J* is a vessel iff 03 is. And as we mentioned 
before, 03** is the same as 53 after a trivial change of coordinates x — > — x on the state space. 
Remarks: 

1. It is easy to check that a^Q^ 1 * (X,t 2 ,t 2 ) and \&(— X,t 2l t 2 ) satisfy the same differential equation. 
Since they are fundamental matrices it is possible iff 

o-^*$- u (\,t 2 ,t° 2 ) = *(-\,t 2> 4)a^* (5.5) 
at u *- u (\,t 2 A) = y*(-X,t 2 ,t° 2 )o-^* (5.6) 

2. This relation between <I> and $ fundamental matrices means that the following relation between 
transfer functions has to be satisfied: 

S(\ > t 2 ) = <T^S:(-\,t 2 )<r 1 (5.7) 

One can easily verify this formula directly, using the vessel conditions and the formulas for S(X, t 2 ), 

s:(-x,t 2 )- 

<rrX*(-A,i 2 M - ar* 1 [D* - B*(-XI + AD^CV;.] V = 

= a^[D - a u C{-XI + A 1 )- 1 B]ax = a^Dax - C{-XI + A^ 1 Ba l = 
= D + C(XI - A!)- x Bai = S(X, t 2 ). 

3. For values of A outside of the spectrum of Ai{t 2 ) (which is independent of t 2 ), we claim that 
S(X,t 2 ) is invertible iff S*(—X,t 2 ) is and from (|5.7|) we conclude that 

S(\,t2)- l = aZ 1 SZ 1 *(-\,t2)a u 



21 



6 System and vessel operations 



As in the classical case we develop basic system operations in our setting, which will be of great impor- 
tance in solving factorization problems for transfer functions. 



6.1 Cascade connection of systems 

Suppose we are given two vessels 

3^' = (A[(t 2 ),F'(t 27 4),B'(t 2 ),C'(t 2 ),D'(t 2 ),D'(t 2 ); 

°i (t a ), a 2 (i 2 ), 7 '(i 2 ), a' u (t 2 ), vL(t2K(t2);H' t2 , £', £'*, £', £'J, 
m" = (A'{(t 2 ),F''(t 2 ,4),B''(t 2 ),C''(t 2 ),D''(t 2 ),D"(t 2 y, 

<(i 2 ), <r 2 '(t 2 ), i'ihWlih), o'L{t 2 )i:{t 2 );H'l 2 , £", £>:, £", £':). 

(6.1) 

and the corresponding systems IT,', IT", defined in (|3.1|) . We want to generate a new system IT by 
feeding in the output of the first system IT' as the input for the second system IT" . To this end we 
assume the output spaces of the first system are the same as the input spaces of the second system (as 
in the case of the classical cascade connection) but also that the corresponding compatibility conditions 
hold: 

o-Uh)=°"(h), ^(t 2 )=a'i(t 2 ), 7 Uh)=l"(t 2 ), £' = £", £'*=£': (6.2) 
Thus we obtain the following system of equations 



IT' 



■£-x'{t u t 2 ) = A'^yx'ih,^) + B'i^y^uit^h) 

x'(h,t 2 ) = F'{t 2 ,tl)x'{ti,tl) + J F'(t 2 ,s)B'(s)a' 2 (s)u(t 1 ,s)ds 

t° 

y'{t u t 2 ) = C'(t 2 )x'(h,t 2 ) + D'(t 2 )«(*i,ta) 



^(ti,t 2 )=^ , (ta)^(t 1 ,* 2 )+fl w (t 2 K.(t 2 )u"(t 1 ,t 2 ) 



IT" 



x"{h,t 2 ) = F"(t 2 , t$)x"(ti,t%) + J F"{t 2 , s)B"(s)a' 2 ^(s)u"(ti, s)ds 

t° 

y(ti,t 2 ) = C"(t 2 )x"(ti,t 2 ) + D"(t 2 )u"{t 1 ,t 2 ) 
Setting u"(t 2 ) = y'(t 2 ), eliminating it and simplifying we get 



IT : { 



a 


' x'ihM) ' 






x"(h,t 2 ) 




' x'(h,t 2 ) ' 






. x"{h,t 2 ) 















x"{ti,t 2 ) _ 


+ 



+ / 

Co 



A[(t 2 ) 
B"(t 2 )a' l (t 2 )C'{t 2 ) A'{(t 2 ) 

p'(t 2 jp\ 

J^F"(t 2 , S )B"( S y 2 \s)C'( S )F(' S ,t 2 )d S F"(t 2 ,4) 

F'(t 2 ,s) 
f; F"(t 2 ,y)B"( y y 2 (y)C'(y)F'(y, t° 2 )dy F"(t 2 ,s) 



B"(t 2 )D'(t 2 ] 
x'(t u 4) 

x"(h,t° 2 ) 

B"{s)D'{s 



o-[u(ti,t 2 ) 



a' 2 (s)u(ti, s)ds 



D"(t 2 )C'(t 2 ) C"(t 2 ) 



x"(tx,4) 



+ D"(t 2 )D'(t 2 )u(t 1 ,t 2 ) 



22 



Thus the corresponding vessel of this system is the following 



353 



A[(t 2 ) 
B"(t 2 )a[(t 2 )C'(t 2 ) A'>(t 2 ) 



B'(h 



ft F"(t 2 , s)B"( S )ai(s)C'(s)F(' S ,4)ds F"(t 2 ,t° 2 ) 



, [ D"{t 2 )C'{t 2 ) C"(t 2 ) 1 ,D"(t 2 )Di(t 2 ),D"(t 2 )D>(t 2 ); 



B"(t 2 )D'(t 2 ) 

a[(t 2 ),a' 2 (t 2 ) n \t 2 ),aUt2)^l(t 2 ) n ':(t 2 );Hl®n' t ^£\£^ £",8'J) 

. (6 - 3) 

Last evaluation suggests that the system IT, with compatibility conditions u'(t\,t 2 ) = y"(ti,t 2 ) indeed 
corresponds to a vessel 

Theorem 6.1 Given two vessels 353', 353", defined in \6.1\) and satisfying compatibility conditions i6.ty) 

<T'l*(t2)=0-'Kt2), 0- 2 *(t 2 )=a' 2 \t 2 ), 7i(t a )=7"(t2), £' = £", £'*=£': 



their cascade connection 353, defined in &6.3\) is a vessel. Transfer functions S"(A, t 2 ), S"'(A, t 2 ), S(X, t 2 ) 
of the vessels 353' , 353" , 353 respectivley satisfy the following relation 

S(X,t 2 ) = S"(X,t 2 )S'(X,t 2 ) 

Proof: We have already seen that 353 is a vessel, provided 353', 353" are. In order to see the formula 
for transfer functions, we feed in the output y'\(t 2 ) of the first system as the input u'^{t 2 ) for the second 
system (recall definitions in section l3^4|) . Then y'\(t 2 ) = S(X,t 2 )u' x (t 2 ) is the input (u'l(t 2 ) — y'\(t 2 )) for 
the second system and 

y'{{t 2 ) = S"(X,t 2 )u'i(t 2 ) = S"(X,t 2 )y' x (t 2 ) = S"(X,t 2 )S'(X,t 2 )u' x (t 2 ). 

We conclude that the transfer function S(X,t 2 ) for the composite system 7E is simply the product of 
the transfer functions of the component systems: 

S(X ) t 2 ) = S'(X,t 2 )S(X ) t 2 ). (6.4) 

□ 

The following theorem is the analogue of theorem 16.11 for differential vessels. It can be either 
deduced from theorem 16.11 by differentiation, or the formulas can be obtained directly by writing the 
systems equations in the differential form and the result established directly by algebraic manipulations 



Theorem 6.2 Suppose that we are given two differential vessels 

2)53' = (Ai(t 2 ), A' 2 (t 2 ), B'(t 2 ),C'{t 2 ), D'(t 2 ),D'(t 2 ); 

aife),^^)^'^),^,^),^^)^:^);^',^',^,^,^) 

and 

2)53" = (A'{(t 2 ),A , ^(t 2 ),B"(t 2 ), C"(t 2 ), D"(t 2 ),D"(t 2 ); 

^fe),^'^)^"^),^^),^^)^"^)^:'^);^"^",^'^",^'). 
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satisfying the compatibility conditions id.ty) , then the following collection 
£2J= ( 

[ D"(t 2 )C'(t 2 ) C"(t 2 ) ],D"(t 2 )D'(t 2 ),D"(t 2 )D'(t 2 ); 



A[(t 2 ) 
B>{{h)o- x {t 2 )C'{t 2 ) A'{(t 2 ) 

B"{t 2 )D'(t 2 ) 



A' 2 {t 2 ) 
B'{{t 2 )o- 2 (t 2 )C'(t 2 ) A>'(t 2 ) 



^(^) ) ^(i 2 ),7'(i2),<(*2),<(t2),7: / (t2)(t 2 );W , ew'',^ J f:, ) £ / ',^) 

is a vessel called cascade connection of the vessels SQJ', . The transfer functions of the corresponding 
systems satisfy the formula Jff.^P - 

6.2 Inversion of systems 

For the classical case [BGKj . if the feed through operator D is invertible, then one can define an inverse 
system having a transfer function equal to the reciprocal of the transfer function of the original system. 
The analogue for t\ invariant overdetermined 2D system is as follows. Suppose that for the vessel in 
the differential form 

®W=(A 1 (t 2 ),A 2 (t 2 ),B(t 2 ),C{t 2 ),D{t 2 ),D{t 2 );a 1 (t 2 ),a 2 ^ 

both D : £ — > £* and D : £ — > £* are invertible. Then we may solve u in terms of y from the last system 
equation 

^x(t u t 2 ) = A 1 (t 2 )x{t 1 ,t 2 ) + B(t 2 ) o-i (<2) u(t u t 2 ) 
: ( -£-x{ti,t 2 )=A 2 (t 2 )x(ti,t 2 )+B(t 2 ) a 2 (t 2 ) u(t u t 2 ) 



by 



y{t u t 2 ) = D(t 2 )u(t u t 2 ) + C(t 2 )x(t!,t 2 ) 



i(t u t 2 ) = -D- 1 (t 2 )C(t 2 )^(ti ! t2) + -D- 1 (t 2 )y(ti,t 2 ) 



and plug it back to get a system _D£ X having the property that (it, a:, y) is a trajectory for DY> if and 
only if (y, x, u) is a trajectory for Z)£ x : 



W[x(h,t 2 ) 
■^-x(ti,t 2 ) 
u(h,t 2 ) = 



= {Ax{t 2 ) - BaD-^t^Cih^xih, t 2 ) + B(t 2 ) aiit^D- 1 ^) y{t u t 2 ) 
= {A 2 (t 2 ) - Ba 2 D- x {t 2 )C{t 2 ))x{t u t 2 ) + B(t 2 ) o 2 {t 2 )D- l {t 2 ) y(h,t 2 ) 
-D-\t 2 )C{t 2 )x{t u t 2 ) + D- X (t 2 )y(t u t 2 ) 



The linkage conditions (|2.10.Link[) means that 



so that this system can be rearranged somewhat to 



irV 2 * =cr 2 D-\ 



DT,' 



^x{h,t 2 ) = (Ai(*a) - B^D-^t^Cit^xihM) + B{t 2 )D- l {t 2 )a ltf {t 2 ) y(t u t 2 ) 
^x{t x ,t 2 ) = (A 2 (t 2 ) - B^D-^t^Cih^xiUM) + B{t 2 )D- l {t 2 )o- 2 *{t 2 ) y(h,t 2 ) 
u{h,t 2 ) = -D- x {t 2 )C{t 2 )x{ti,t 2 ) + D- 1 (t 2 )y(t 1 ,t 2 ) 



This suggests the following theorem 
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Theorem 6.3 The following collection 
W x =(A 1 (t 2 ) x , A x (t 2 ), 5 x (i 2 ), C x (t 2 ), ^ x (t 2 ),5 x (t 2 );a 1 ,(i 2 ), ( T 2 ,(t 2 )7 >t (i 2 ),a 1 fe), ( 7 2 (t 2 ),7(t 2 ); 

where _ _ 

4i(t 2 ) x = Mt 2 ) - BaxD-\t 2 )C{t 2 ), A x (i 2 ) = A a (i a ) - Ba 2 D-\t 2 )C{t 2 ), 
B x (t 2 ) = B{t 2 )D-\t 2 ), C x (t 2 ) = -D-\t 2 )C{t 2 ), 
D*{t 2 )=D-\t 2 ),D*{t 2 ) = D- 1 {t 2 ) 

is a vessel with transfer function S x (X,t 2 ) equal to the inverse of the transfer function S(X,t 2 ) of the 
vessel 23. 

Proof: We have to show that all the vessel conditions hold. Let us omitc the t 2 dependence of all the 
operators 

• Lax equation: gpA* = A* A% (t 2 ) - A| (t 2 )A* (t 2 ). 

W 2 A i =M A ^~ Bo-iD-iC] = £A, ^[BaiD-iC] = ±A X £[Bax]D^C - Ba^D^C] 
On the other hand, 

A x A x (t 2 ) - A x {t 2 )A x (t 2 ) = {Ax - Bo-iD~ x C) (A 2 - Ba 2 D- l C) - [A 2 - Ba 2 D~ l C){Ax - Bo\D~ l 
= A X A 2 - A 2 Ax + 

+[-A x Ba 2 + A 2 Bo- 1 - Bo 2 D- 1 CBo 1 ]D- 1 C+ 

+Ba 1 [-D- 1 CA 2 + a^ 1 cr 2 D- 1 CA 1 + D^CB^D^C] 

Using now linkage conditions (|2.10.Link| for the original vessel 2)23 the result follows. 

• Input vessel condition: ^■(i? x c r i*) - A^B x o x * + A x B x a 2 * + -B x 7* = 0. We first simplify the 
expression with derivative: 




and the other elements are 

-A x B x a 1 ^+A x B x a 2 ^+B x ^ = 

= -[M- Ba 2 D-^C]BD-\r u + [Ax- BaxD^^BD^a^ + BD~ x ^ m = 
= \-A 2 Bo x + A x Bo- 2 - BD- 1 o- 2 *CBcj 1 - BD^o^CBo^D- 1 + BD^ 1 ^ 

arid using now the linkage condition we obtain that their sum is zero. 

• Output vessel condition: ^(ciC x ) + eriC x A 2 + o 2 C x A x + jC x = is similar to the input 
vessel condition. 

• Linkage conditions: 

oxD x = D x ax*, o 2 D x = D x a 2 ^ 
D x ^=a 2 C x B x a 1 ,-<7 1 C x B x a 2 ,+a 1 ^ ; D x +^D X . 

This is an immediate result of the linkage condition for the original vessel, rearranging the elements 
and multiplying by inverse of D. 
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□ 

The analogue of the last theorem for vessels in the integral form can be obtained by integrating 
the corresponding formulas of differential vessels. The result is as follows 

Theorem 6.4 Given the integral vessel 35J 

3^={A 1 (t 2 ),F(t 2 ,t 2 ),B(t 2 ),C{t 2 ),D{t 2 ),D{t 2 )] 

(£2), o- 2 (t 2 ), 7(^2), o-i* (t 2 ), 02* (£2)7* (^2); Ti t2 , £, £*,£, £*) 

the following collection is a vessel (called inverse) 

3% x = (A?(t 2 ),F*(t 2 ,t° 2 ), B(t 2 )b-\t 2 ),-D-\t 2 )C{t 2 ),D- l {t 2 ),b- l {t 2 )- 

o-i*{t 2 ),a 2 ^(t 2 )^^(t 2 ),a 1 (t 2 ) 1 a 2 (t 2 )^(t 2 )\H 1 £*,£,£*,£), 

where 

Ai(i 2 ) x = Ai(t 2 ) - Sa 1J D- 1 (t 2 )C(t 2 ), 

F>< (t a , i°) = F(ta, t§) - /| F(t a , 4)B(4)a 2 (^-^^(4)^(4 ^ dt\-...- 

F(tltl)---F(^-\q)B(q)a 2 (q)D-\^ ■■■dt 1 2 -.... 

Proof: Notice that the formula for F^,^) ^ s j us ^ Peano-Baker formula for F(t 2 ,t 2 ) generalized from 
the differential equation (holding for the differential vessel) 

J-F x (t 2 , = A x (t 2 )F x (t 2 ,t°) = -L4 2 (t 2 ) - Ba.D-Ht^CihKF* (t 2 ,t° 2 ) 
at 2 

□ 

6.3 Projection, compression and cascade decomposition of systems 

Following the construction of cascade connection, it is natural to ask whether the reverse construction 
exists. One of the main ingredients of this construction is that the state space Ht 2 is decomposed into 
two subspaces H = TL' t2 © Tt" 2 , which are invariant for the following operators: 

Ai(i a )7^ a C7^ a , F{t 2 ,t° 2 )n' t0 =n' t2 , Vt 2 ,t° 2 
A? (t 2 )H>> 2 C HI , F>< (t 2) 4)H>;» = H>> 2 , Vt 2 ,4 

In the differential case this means that H' t2 is invariant under A\{t 2 ) , A 2 (t 2 ) for all t 2 and that H' t ' 2 is 
invariant under A* (t 2 ), A 2 (t 2 ) for all t 2 . 

Definition 6.1 Given a vessel 35J subspaces Gt 2 £ Ht 2 are called invariant if for all t 2 holds 

A x {t 2 )g t2 cg t2 , F(t 2 ,t° 2 )g t o =g t2 . 

Subspaces C? t x C Tlt 2 are called co-invariant if their complements Ttt 2 Qgf 2 are invariant and if these 
subspaces are invariant for the operators A* (t 2 ) and F x (t 2 ,t 2 ) (defined in theorem \6.J$ , i.e., if for all 
t 2 holds 

A*(t 2 )g*cg*, F*(t 2 ,t° 2 )g*o=g? 2 . 
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The classical condition for a cascade decomposition and a factorization of the transfer function [BGK, 
BFKD, SJ uses these two notions of invariance. We present an analogue of the corresponding theorems. 
Assume that we are given an overdetermined 2D system, t\ invariant (|3.1[) with the vessel 

303 = (Ai(t 2 ), F(t 2 , 4), B(t 2 ), C{t 2 ) 7 D(t 2 ), £>{t 2 ); CTl (t 2 ), a 2 (t 2 ), j{t 2 ), tru(t 2 ), ^.(tafr.te); 

Suppose also that we are given subspaces (? t2 C 7i t2 that are invariant. Then it is possible to define a 
projection of the vessel 303 onto the invariant subspaces Gt 2 as follows 

Definition 6.2 Projection of the vessel 303 on the invariant subspaces Gt 2 * s a collection 

303' = (A[(t 2 ), F'(t 2 , B'(t 2 ), C'(t 2 ), D(t 2 ), D(t 2 ); <n(t 2 ), a 2 (t 2 ), 7 (i a ), ^(t 2 ), <(t 2 )7l^2); ^i 2) ^ , 

where denoting by Pg t2 - projection on Q t2 , the operators in 303' are 

F'(t 2 ,t 2 ) = F(i 2 ,f 2 )P G ( t o), 
A' 1 (t 2 ) = F'(t 2) ^)A 1 (*0)F'(tO,t 2 ), 
B'(t 2 ) = Pg t2 B(t 2 ), 
C'(t 2 ) = C(t 2 )Pg t2 , 

and <T^(t 2 ), cr 2 *(t 2 ),7*(i 2 ) are taken so that the linkage conditions t2.10.Link^) are satisfied 

a' u D = Da 1 , a' 2 *D = Da 2 , D 7 = a'^C'B'a^ - a[*C'B'a 2 + a' 1)f ^-D + jiD. 

dt 2 

Let C/ t2 be co-invariant subspaces. Then we define compression of the vessel 303 onto the the co-invariant 
subspaces Q x 2 as follows 

Definition 6.3 Compression of the vessel 303 on the co-invariant subspaces Q x 2 is a collection 

303* = (Af^),^-^,^),^-^),^^^),^^),^^); _ 
^'(*2W(i 2 ),7"(t 2 ), «ri.(t a ), cr a ,(t a )7.(*a); 0£, £.)• 

where A± (t 2 ), F x (t 2 , t 2 ) are defined as in theorem \6.4\ and denoting by P x x - projection on Q x 2 

B*(t2) = P^Bfe), C"(t 2 ) = C(t 2 )P e x t2 , 
</ie operators <r"(t 2 ), a 2 (i 2 ),7"(i 2 ) are taken so that the linkage conditions \2.10.LinM) are satisfied 

truD = Do / S, o- 2 ,D = Da' 2 ', 5 7 " = a 2 „C x B x of - a u C x B x a'i + tr u -^-D + j,D 

dt 2 

Lemma 6.5 303',303 x are vessels. 

Proof: We will show that 303' is a vessel. For 303" the proof is essentially the same. 

First we show that F'(s, t) is an evolution group. F'(s, s) = F(s, s)P g ^ s - ) = Ig(t°), 

F'(s, t)F'(t, y) = F(s, t)P g(t) F(t, y)P g{y) = F(s, t)F(t, y)P g(y) = F(s, y)P g(y) = F'(t, y). 
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Then we have to show that all vessel conditions are satisfied. 

A[(t 2 )F'(h,t 2 ) = F\h,tl)A x {tl)F\t^)F\h,t%) = F{h^)A x (? 2 )P g{tl) =F\h^)A[(? 2 ), 
which means that the Lax equation holds. In order to check the input vessel condition notice that 

F'(t° 2 ,t 2 )B'(t 2 ) = P g{t o } F(t 2 ,t 2 )P gt2 B(t 2 ) = P G{t o 2) F(t Q 2 ,t 2 )B(t 2 ), 
and in the same manner 

F'(4, t 2 )A[(t 2 )B'(t 2 ) =F'(t° 2) t 2 )F'(t 2 , t^A^F'itl t 2 )P St2 B(t 2 ) = 

= Al(t^)Pg(_ t OjF(t^, t 2 )B(t 2 ) = Pg( t O)F(t 2 ,t 2 )Al(t 2 )B(t 2 ). 

So, the input vessel condition is 

^(F'itlhWih)^)) + F'(t° 2 ,t 2 )A' 1 (t 2 )B'(t 2 )a 2 {t 2 ) + F'(t a 2 ,t 2 )B'(t 2 ) Z (t 2 ) = 

= PQ&)[&(F(4,t2)Bfa)<ri(t2)) + F{t° 2 ,t 2 )Ai(t 2 )B{t 2 )a 2 (t 2 ) + F(t° 2 , f 2 )B(t 2 h(f 2 )] = P G (t°)0 = 0- 

The output vessel condition is a consequence of other conditions by considering the differential equation 

J-5'(A,t 2 ) - o^{t 2 ){Xa 2 *{t 2 ) +iM))S'(\t 2 ) ~ S'{\t 2 )°i 1 {t 2 ){\a 2 {t 2 )+ 1 {t 2 )) 
ot 2 

for the projected transfer function 

S'(X,t 2 ) =D(t 2 ) + C'(t 2 )(XI-A' 1 (t 2 ))- 1 B'(t 2 )<j 1 = 

= ^(t2) + C(t 2 )F(t 2 ,^)P t o(A/-A 1 (tO))- 1 P t oF(t°,t 2 ) J B(i 2 ) ( T 1 . 1 ■ ' 

Let us first evaluate the derivative 

£S'(A, t 2 ) = i- 2 D + ^[C(t 2 )F(t 2 ,t° 2 )P t o](XI - A 1 (4))- 1 PqF(l%,t 2 )B(t2)a 1 + 

+C(t 2 )F(^,i§)P t o(A/-A 1 (tO))- 1 P t o^[F(i«,i 2 )i3(t 2 ) ( 7 1 ] = 
= + a^^C^F^Mt^o + j'C(t 2 )F(t 2 ,t° 2 )P t o}(XI - A 1 (t° 2 )r 1 P t oF(t° 2 ,h)B(t 2 )a 1 - 

-C(t 2 )F(t 2 ,t° 2 )P t o(XI - A 1 (t0))- 1 P t o[A 1 (^)P(t°,t 2 )i?(t 2 )a 2 +F(tO,< 2 )i?(t 2 ) 7 ] = 
= ^-D + a^[a 2 *X + 7']C(i 2 )P(t 2 ,t°)P t o](A7 - A 1 (t° 2 ))-^P t oF(tl t 2 )B{t 2 )ax- 

-C(t 2 )i!'(t 2) ^)P^(A/-A 1 (^))- 1 P^ 1 (^)i!'(^,t 2 )B(t 2 )[Aaa + 7]- 

-^a 2 ^C{t 2 )F(t 2 ,t%)P t oF{tl,t 2 )B^ 
= £- 2 D + a^[a 2 *X + i](-D + S'(X, t 2 )) - (-D + S'(X, t 2 ))a x l [Xa 2 + 7]- 

-Vu 1 (T2Mt2)F(t2,4)P t aF(tlh)B(t 2 )cT 1 + C (t 2 )F(t 2 , t° 2 ) P t oF(t° 2 , t 2 ) B (t 2 )a 2 = 
= ^D + erf* V 2 *A + i]{-D + S'(X, t 2 )) - (-D + S'(X, t 2 ))a^[Xa 2 + 7]- 

-a 1 -V 2 ,C'(t 2 )P t oP'(t 2 ) C 7 1 + C'(t 2 )P t oB'(t 2 )a 2 . 

If we plug this expression into (|6.5p . we shall obtain 

£- 2 D - [a 2 *X + i]D + Da^[Xa 2 + 7]- 

-a i ;V 2 *C ,, (t 2 )P t aP'(f 2 )(7i + C'(t 2 )P t oB'(t 2 )a 2 = 

□ 
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It is also possible to perform a compression of the vessel on a semi-invariant subspace. Let us 
first define it. Suppose that Gt 2 is an invariant subspace of the vessel 303. Suppose also that there is a 
co-invariant within Gt 2 subspace G[ 2 C Qt 2 , that is the subspace Gt 2 —G't 2 is co-invariant. In this case we 
call Q' t2 a semi-invariant subspace of 303. Then performing projection of the vessel 303 on the invariant 
subspace Gt 2 > we shall obtain a new vessel 303 . Performing further compression on the co-invariant 
subspace G' t2 we shall obtain the desired vessel 303". 

Remarks: 1. Notice that invertibility of B(t 2 ) is not essential but suffices in order to determine all 
the relevant data. For example, in the case of projection on an invariant space Gt 2 one obtains 

</ u = Da x B-\ <4 = Da 2 D-\ ^ = BjB' 1 - \o'^C'B'o x - a' u C'B'a 2 + a^^-BjB' 1 

dt 2 

and consequently they are uniquely determined. 

2. It is enough to have a subspace Q t o C H t o that is invariant under Ai(t 2 ) with a complementary 
subspace £7 t * C 7i t o that is invariant under A* (t 2 ) and such that 

F(t 2 ,4)g t o, F(t 2 ,t° 2 )g* 

2 L 2 

are complementary for all t 2 . 

3. For the conservative case (to be studied later) A* (t 2 ) = A\{t 2 ), so the existence of a complementary 
invariant subspace is automatic, just like in the classical case. 

As a result of all these consideration we deduce a theorem of cascade decomposition of vessels. 
Suppose that we are given an invariant subspace Gt 2 , which is at the same time co-invariant. Then it 
is possible to produce a projection on Gt 2 and compression on its complement. The point is that the 
obtained in such a way vessel can be cascade connected to give the initial one. This is precisely the 
content of the next theorem 

Theorem 6.6 Suppose that we are given a vessel 303 of the form 16.3]) with invariant subspace Qt 2 , 
and co-invariant subspace G* 2 , then the projection on Gt 2 produces a vessel 303' and compression to Gf 2 
produces 303 x . Moreover, it is possible to cascadly connect these two vessels and to obtain the original 
one 303. 



7 Kalman decomposition of Vessels 



The notions of approximate controllability and observability allows building of minimal systems for 
which there is a very good classification theory. It turns out that there are possible other parts of a 
system (vessel) that are non-approximatcly controllable or are not observable. Let us denote the system 
(|33]1 (for fixed a x (t 2 ), a 2 (t 2 ), 7(i 2 ), 0i*(t 2 ), <?2*(t 2 ), 7*(*2)) 



IT, : < 



by 



£x(h,t 2 ) = A x (t 2 ) x(t!,i 2 ) + B{t 2 )a 1 (t 2 ) u{tx,t 2 ) 

x{h,t 2 ) =F{t 2 ,t 2 )x(t 1 ,t° 2 ) + J F(t 2 ,s)B{s)a 2 (s)u(t 1 ,s)ds 

t° 

K y(h,t 2 ) = C(t 2 ) zfa.taJ + DSaMti.ta). 

IT = [Ai(£ 2 ), -F(£2i t 2 ), B{t 2 ), C(t 2 ), B(t 2 )]. 
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Theorem 7.1 Let IT be a system defined in \3. 1\) with inner state space Ttt 2 for each t 2 . Then there 
exists an orthogonal (for each t 2 ) decomposition of the state space 



Hi 



H c t °< 



>H c t 



n 



1-2 



so that with respect to this decomposition the system has the following decomposition of its operators 

A x (t a ) = 



Aia(t 2 ) 

A?{t2) 






B{h) ■ 
C(t 2 ) = 
F(t 2) t° 2 ) = 



Af(t 2 




B co 

B co 



C co (t 2 ) C 5 °(t 2 ) 
" F cB (t 2 ,f 2 ) ~ 



A 13 (t 2 ) 

A 23 {t 2 ) 

Af°(h) 




A M (i 2 ) 
A 2i (t 2 ) 
A M {t 2 ) 
A\°{t 2 ) 



F c °(t 2l t° 2 ) 

F™(t 2 ,4) 
F eo (t 2 ,t° 2 ) 



Af{t 2 ) 




A 12 {t 2 ) 
A?{t 2 ) 



where the subsystem, defined by 
IT C = 

is approximately controllable, the system 

IT° = [ 
is observable, and the system 



Al°{t 2 ) A 24 (t 2 ) 
Af{t 2 ) 



F c5 (t 2 ,t° 2 
F c °(t 2 ,t° 2 



F co (t 2 ,t° 2 ) 
F eo (t 2 ,t° 2 ) 





B co 


1 









[ C c °{t 2 ) } ,D] 





■ B co ■ 


1 






,[ C co (t 2 ) C s °(t 2 ) ],D] 
[Af ,F co {t 2 ,t° 2 ),B co (t 2 ) , C co (t 2 ) , D] 



is minimal (i.e., approximately controllable and observable). Moreover, the transfer functions of all the 
systems are equal: 

s(\,t 2 ) = s c {\t 2 ) = s°(\,t 2 ) = s co (\,t 2 ). 



Proof: Denote for each t 2 



g c (t 2 ) = spm{im A{(t 2 )B(t 2 ) I j = 0,1,2,...}, g s (t 2 ) 

g o (t 2 ) =span{Im^'(t 2 )C*(t 2 ) | J = 0,1,2,...}, G B (t 2 ) 



n t2 eg c (t 2 ) 
n t2 eQ {t 2 ) 



Now write each operator in the system IT as an operator with respect to Q c {t 2 ) and Gs{t 2 ). Since c (t 2 ) 
is Ai(t 2 ) invariant, the system IT = [Ai(t 2 ), F(t 2 , t 2 ), B(t 2 ), C{t 2 ), D{t 2 )\ is of the following form 



,[ C c (t 2 ) C 5 (t 2 ) ],D(t 2 )} 



' A\{t 2 ) A 12 (t 2 ) ' 




' F°(t 2 ,t° 2 ) - 




" B c (t 2 ) ' 


A\{t 2 ) 


5 


. F%t 2 ,t° 2 ) _ 


5 
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and clearly 

C c (t 2 )A c (t 2 )B c (t 2 ) = C{t 2 )A{t 2 )B(t 2 ). 
Thus the original system IE has the same transfer function as 

7£ c = [Al(t 2 ),F c (t 2 ,t° 2 ),B c (t 2 ),C c (t 2 ),D(t 2 )] 

and this system is approximately controllable. The same process works on the given system (|3.1|) with 
observability (Q 5 rather than Q is invariant under A\{t 2 )) to give a "smaller" observable system 

IE" = [Al(t 2 ),F%t 2 ,t°),B°(t 2 ),C°(t 2 ),D(t 2 )]. 

If one combines these two processes, one gets (by first decomposing the system IE into controllable and 
uncontrollable parts and then decomposing these systems into observable and unobservable parts) the 
desired decomposition. □ 

8 Analytic functions as transfer functions 

The aim of this section is to show that any function in our class X can be realized, i.e., presented as a 
transfer function of a certain vessel. 

8.1 Realization theorem for arbitrary analytic functions in X 

So, suppose that we are given a function S(X,t 2 ) £ X. Our first aim is to realize this function, i.e. to 
show that this class is realizable. In order to do it, we realize S(X,t 2 ) for a fixed t 2 [H] as 

S(X,t° 2 ) = D + C (XI- A 1 )- 1 B a a l {tl). 

Then the following theorem holds 

Theorem 8.1 Suppose that S(X,t 2 ) £ X. Then there exists vessel 2)03 in the differential form 

®y={A 1 (t 2 ),A 2 (t 2 ),B(t 2 ),C\t 2 ),D(t 2 ),D;a 1 ,a 2 ^,a 1 *,a 2 *,^;H,£, £*,£,£,), 
with this transfer function and for which 

C(ta)= j> $*(X,t 2 ,t 2 )C (XI - A^dX (8.1) 

SpecAi 

B(t 2 )= I {XI-A l )- 1 B a 1 {t 2 W{-XMA)^ (8.2) 

SpecAi 

D{t 2 ) =S(oo,t 2 ) (8.3) 

and 

S(A, t 2 ) = D(t 2 ) + C(t 2 )(XI - Ax)- x B{f 2 )(Tx. 

Before we consider the proof of this theorem, it is important to note that for the functions C(t 2 ), B(t 2 ) 
the following lemma holds 
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Lemma 8.2 C(t 2 ), B(t 2 ) satisfy the following differential equations with the spectral operator parameter 
Ax: 

^ [B(ta)cri(ia)] = A 1 B(t 2 )a 2 (t 2 ) + B(t 2 ) 7 (i a ) 
Proof: For C(t 2 ) one obtains that 

o-iC(t 2 y =<j 17 ^ § $*(\,t 2 ,4)c (\i-A 1 )- 1 d\ = 

SpecAi 

= <n § <Ji 1 (a 2 \ + ^)<P*(\,t2,t 2 ) )Co{\I - A 1 )- 1 d\ = 

SpecAi 

= § a 2 ^,{X,t 2 ,t° 2 )C (XI - A x + A 1 ){XI - A 1 y 1 d\+ 

SpecAi 

+ § 7,$»(A,t 2 ,t§)Co(A/ -A 1 )- 1 d\ = 

SpecAi 

= a 2 § ^ if {\t 2 ,t° 2 )C {XI-A 1 )- l dXA l + ^ § ^(X,t 2 ,4)C (XI -A 1 )- 1 dX = 

SpecAi SpecAi 

= a 2 C(t 2 )Ai + j*C{t 2 ) 

and the same proof works for B(t 2 ). □ 
Proof of theorem I8.lt Let us define a vessel: 

2)93 = (A 1} A 2 = 0,B(t 2 ), C(t 2 ),D(t 2 ),D = a^Dfe)^ 1 ; a 1} a 2 , 7, a u , cr 2 *, 7*; W, £, £*, £, £*), 

and show that all vessel conditions hold. Lax equation holds 

4-Ai = = A X A 2 - A 2 A X = A X Q - OAl 
dt 2 

The first (|2.5.0vcrDp and the second (|2.9.0vcrD[) vessel conditions are exactly the contents of lemma 
ET21 

Consider next the expression C(t 2 ){XI - A 1 )- 1 B(t 2 )ai(t 2 ). Using lemma (JO) we obtain that 

i- 2 {C(t 2 )(XI - A^Bit^o-iih)) = 
o-u 1 (* 2 *C(t 2 )A 1 + ^(t 2 )C(t 2 ))(XI - A 1 )- l B{t 2 )a 1 {t 2 )+ 

+C(t 2 ){XI - A 1 )- 1 [A 1 B{t 2 )a 2 {t 2 ) + B(t 2 )j(t 2 )] = 
= o-^(a 2 *X + j*(t 2 ))C(t 2 )(XI - A^Bit^ih) - C(t 2 ){XI - A 1 )- 1 J B(t 2 )(a 2 (t 2 )A + 7 (t 2 ))- 

We can use here the fundamental matrices $(A, t 2 , t®), $*(A, t 2) t®) and with the help of variation of 
coefficients, we obtain that 

C{t 2 ){XI - A 1 )- x B{t 2 )a l {t 2 ) = ^{X,t 2 ,t° 2 )K(X,t 2 )^- 1 (X,t 2 ,t° 2 ), 

where 

^K(X,t 2 ) = 

$-HA,i 2 ,<°)(-a 1 -> 2 ,C(t 2 )S(i 2 )a 1 (t 2 ) + C'(t 2 )S(i 2 )a 1 (t 2 )^ 1 (i 2 )a 2 (t 2 ))$(A^ 2 ,^ 
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with the initial value i^(A,t§) = C^XA/ — A\) 1 B(t 2 ) )tJi(t 2 \). Since the fundamental matrices are 
entire in A functions, we also obtain that 

K(X,t 2 )-K(X,t°) = [* ^-K(X,y)dy 
Jt° "V 

is an entire function of A. Thus 

S(X,t 2 )-C(t 2 )(XI - Ai^Bfafafa) = ^(X,t 2 ,4)(S(X,t° 2 )-K(X,t 2 ))^- 1 (X,h,t° 2 ) = 
= $,(A,t 2 ,ti)( J D(^)+C(tO)(AI-A 1 )- 1 J B(tO)a 1 (^)- 

-C(t° 2 )(XI-A 1 )-^B(t 2 )a 1 (t° 2 )- f%%K(\,y)dy)Z- l (\,t 2 ,4) = 
= *,(A,t 2 , *§)(£>(*§) + ^K(X,y)dy)^- 1 (X,t 2 ,t° 2 ), 

which is an entire in A function too. On the other hand, when A tends to oo, 5(00,42) = D(t 2 ) and 
lim C(t 2 )(XI — Ai)~ x B(t 2 )a\(t 2 ) = 0, which means that their difference is bounded and consequently, 

A — >oc 

by Liouville's theorem for operator valued functions, applied to constant vectors is constant. Finally, 

S(X,t 2 ) =D(t 2 ) + C(t 2 ){XI-A 1 )- 1 B(t 2 )a 1 (t 2 ). (8.4) 

Once we have established all these formulas, it remains to show that the linkage conditions (|2.10.Link|) 

(Ti*D = Dai, <J 2 *D — Da 2 

Lh) = a 2tf CBa x - a^CBa 2 - a u + 

are satisfied. In order to do this, we use the differential equation for S(X,t 2 ) (|3.1ip 

J-5(A,t 2 ) -a 1 - 1 (t2)( ( T2*(t 2 )A + 7>t (t 2 ))5(A,t 2 )-5(A,t 2 )ar 1 ^2)(a2(t2)A + 7(t2)). 
ot 2 

and substitute here the realization formula (|8.4[) . Then 

±D(t 2 ) + i-C{t 2 ){XI - Ai)~ l B(t 2 )ax{t 2 ) + C(t 2 )(XI - A^^B^)^)) = 
= o^{t 2 ){a 2 *{t 2 )X + j*(t 2 ))(D(t 2 ) + C(t 2 )(XI - A 1 )~ 1 B(t 2 )a 1 (t 2 ))- 
-(D(t2) + C(t2)(A/-A 1 )- 1 J B(t 2 )a 1 (t2))^r 1 (i2)(a2(t2)A + 7(i2)). 

Considering the linear in A part, for big values of A we immediately obtain that 

and plugging this back into the differential equation, and tending A to infinity, we obtain (defining 
D = a^Da^ 1 ) that 

Dj = a 2 ^CBa\ — a\^CBa 2 — a±* — — D + 7*£>, 

dt 2 

which finishes the proof. □ 

Notice that there are no assumptions on the dimensions of £, £*. In the next section we consider 
the finite dimensional case dim7i t2 < oo. 
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8.2 Realization theorem for matrix functions in X 

In this section we want to further investigate the formulas for C(t 2 ),B(t 2 ) arising in the realization of 
S(X,t 2 ) €E K for the matrix case. So, S(X,t 2 ) maps solutions of the input ODE (|3 . 7[) with the spectral 
parameter A to solutions of the output ODE (|3.8p with the same spectral parameter. By Jordan theorem 
any constant matrix is similar to its unique Jordan form, and consequently, we are going to concentrate 
on this special case. More explicit realization, based on the theorem 18. II is achieved in the next theorem. 
So, suppose that there is only one eigenvalue for Ai of multiplicity n, 

Theorem 8.3 Suppose that S(X,t 2 ) has one pole z of maximal order n. Then there exists a chain of 
functions {co*, ci*, . . . , c„*} and {bo, bi, ... , b n } such that 



S(X,t2) = D + [ c * cu ... €(„_!)*] (XI - Jordan(z)) 



1. where Co* is a solution of the output differential equation k3. 8\) with the spectral parameter z, and 
Ci* is a solution of the differential equation 

d 

za 2 *y - <ri*— y + 7*2/ = c 2 *Ci_i*. 
ot 2 



"0 

b\ 



2. and bo is a solution of the adjoint output differential equation \5.J$ with the spectral parameter 
—z*, and bi is a solution of the differential equation 

* 9 , 

-z a 2 y - o- 1 —y + = a 2 bi-i. 
ot 2 



Proof: Remember that S(X,t 2 ) is of the form (|3.10|) 

S(X,t 2 ) = $,(A,^,i 2 )^(A,^)$- 1 (A,t°,t 2 ), 

where <&(A, t 2 , i 2 ), <£>*(A, t 2 , t 2 ) are fundamental matrices for the input and the output ODEs, respectively, 
and S(X,t 2 ) is a rational in A matrix, for which we can apply the realization theorem: 

S(X, 4) = D + C(XI - Jordan^])' 1 B 



Then we have 



S(X,tlt 2 = ^,(X,tlt 2 )S(X ,t 2 )^- 1 (X,t 2 ,t 2 ) 

1 1 1 

A-z (A-2) 2 (A-z) 3 

n 1 1 

U A=I (A-z) 2 

y^- 

X — z 



= $*{X,4,t 2 )[D + C 



(A-z) — 



1 

A-2 



B]^(X,t%M)- 
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Writing now the fundamental matrices in the Taylour series we obtain that 

oo oo 

oo oo 

$- 1 (A,i°,< 2 )-<i>- 1 (z,iO,i 2 )+ £ $ fc ,(tO ;i2)(A _ z) fe = £ <& fe ,(A-z) fc 

fe'=l fc'=0 

Inserting these expressions for the fundamental matrices, we are able to calculate the Lourent coefficients 
explicitly: 

1. The coefficient of (x-z) n * s 

where we denoted Co* = <E>o*C e i, &o = o'j~ 1< I ) o-Be Tl . 

2. The coefficient of tt — \„_ t is 

(A-z)" 

^o.Ceie^xfl^o + $o*Ce 2 e*B$ + ^Q*Ce x e* n B^i + $> u Ceie n B$o = 
= Co*[e*-i-B$ + e*S$ x ] + [<J> *Ce 2 + ^i*Ce 1 ]b* cr 1 = 

— Co*6*CTl + Cl*&g(7l 

Notice that direct calculations show that ci*,&i are companion solutions of co*,&cb respectively. 

3. The coefficient of tt — \„_ 2 is 

(A-z)" 

*o*Cei(e;_ 2 B$ + <_i-B$i + e„S$ 2 ) + ($ *Ce 2 + ^CeiXe^-B^o + <#$i)+ 
($ *Ce 3 + $i»Ce 2 + $ 2 „Cei)e nJ B$ = 

= C *6 2 CTl + Ci*^^! + C 2 *&o<7i, 

where we denote fc^i = e£_ 2 -B<I>o + e^_x-B$i + e n B§ 2 and c 2 * = $ *Ce 3 + <I>i*Ce 2 + $ 2 *Cei. 
Again, as in the previous case, it is a matter of simple calculations to show that c 2 *,6 2 are 
companion solutions to ci*,&i, respectively. 

4. For all the other coefficients, by induction, we obtain the desired result. 

Notice that $ fc *(t§,t 2 ) = gfir$*(A,t£,t2)U=», $k{4,h) = -£ J <$>~ 1 {\ 7 4 7 t 2 )\\= z and this is used to show 
the connection between Ck* and C( k+1 ^ (pk and bk+i). D 

It means that one obtains a vessel with the transfer function S(X,t 2 ) by means of the following 
definitions: 

Ai = Jordan(zj), A 2 (t 2 ) = 0, C(t 2 ) = Col{c«(t 2 )}, B{t 2 ) = Row {h(t 2 )} 

It is a matter of simple calculations to show that all the vessel conditions are satisfied. 

As for the general case, suppose that A\ — [/fit/ -1 , where f2 is the Jordan block form of A\. 
Then from realization theorem 18. 11 we obtain that 

S(X,t 2 ) = D(t 2 ) + C{t 2 ){\I - A 1 )- 1 B(t 2 )a 1 = D(t 2 ) + C(t 2 )V-\\I - Sl)- 1 VB(t 2 )<r 1 
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and one can consider the case where A\ is of the Jordan block form. Suppose that we are given chains 
of the appropriate sizes of companion solutions Cq„, . . . , c l ^ n and 6^, . . . , &L._i\„ for the output and 
the adjoint input ODE with the spectral parameter A;. Suppose that it is given for each eigenvalue 
Xi of f2. Then the final vessel is obtained by defining C(t 2 ) as a block matrix, where the matrices 
c o* ■ ■ ■ c \ ni -i)* are 011 ^ ne diagonal. The same construction works for B(t 2 ) 1 so that the final result and 
vessel conditions are easily verified. 

The construction of the system parameters from the residues of the given rational matrix function 
(which one assumes has only simple poles), is known as the Gilbert realization (see [K], page 349). 

It is also appropriate to emphasize here that we have built in theorem 18.31 right and left pole 
pairs for the matrix function S(X, t 2 )- Let us recall first the definitions (from |BGR| ). From the formula 

s(\,t 2 ) = <f>4\,t 2 ,T 2 )s(\,t° 2 )<s>- 1 (\,t 2 ,t° 2 ) 

it follows that the poles of the matrix S(X,t 2 ) are independent of t 2 . So, suppose that z is a pole of 
S(X,t 2 ), which has the following Lourent expansion around z: 

00 

S(X,t 2 )= ]T(A-z)^(i 2 ). 

Then yo{t 2 ), ■ ■ ■ ,y r -i(t 2 ) is called a right pole chain for S(X,t 2 ) at z, if there exist additional vectors 
y r (h), ■ ■ ■ ,Vr+q-i{h) (9 being the order of z as a pole of S^ 1 (X,t 2 )) such that S~ 1 (X 7 t 2 )y(X,t 2 ) is 
analytic at z with zero of order r at z, where 

r q -l 

y(X,h)= ^(A-z)V 
3=0 

The natural number r corresponds to a size of Jordan block (with z on the diagonal) to be constructed 
shortly. 

A canonical right pole pair (A\, B(t 2 )) is a collection of chains 

B{t 2 ) = y W(h), . . .d-ii^VoHh), ■ ■ .,y^ 1 (*2)i ■ • ■ ;y { p \t2), ■ ■ -d-iih); 

and corresponding Jordan blocks 

Ai = Ji © J 2 © ■ • ■ © J p . 

with the property above. Moreover, any pair similar to (Ai, B(t 2 )), i.e., for an invertible matrix M(t 2 ), of 
the form (Mft^AiM -1 ^), M(t 2 )B(t 2 )) is called a right pole pair. Analogously, one defines a canonical 
left pole pair (C{t 2 ),Ai), when y{X,t 2 )S~ l (X,t 2 ) is demanded to be analytic at z with zero of order r 
at z. 

In theorem 18.31 we have constructed such pairs for each t 2 , but even more, the left pole pair 
(C(t 2 ), Ai) has an additional property, that C(t 2 ) satisfies the differential equation in lemma [8?2l with 
the spectral matrix parameter A\. This in turn means that the left pole chains consist of companion 
solutions of the output ODE ()3.8() for poles of S(X, t 2 ) as described in theorem 18. 31 A similar conclusion 
is true for {Ai,B(t 2 )). 

Finally, a triple [C{z,t 2 ), A\{z), B{z 1 t 2 )a\{t 2 )) is called a pole triple at z of a rational matrix 
function S(X,t 2 ) if (C(z, t 2 ), A\{z)), {A\ (z), B(z, t 2 )a\ (t 2 j) are the left and the right pole pairs at z, 
respectively. In this case 

S(X, t 2 ) - C(z, t 2 )(zl - A^z^Biz, t 2 )o- 1 {t 2 ) 
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is an analytic function at z. A (full) pole triple of S(X,t 2 ) is just a direct sum of all local pole triples, 
which respects the order of their appearance: 

C(ta) = C{z u t 2 ) © C(z 2 ,t 2 ) © ■ ■ • © C(z p , t 2 ), 
Ai = diag{Ai(zi),^i(z 2 ), ■ • .,Ai(z p )}, 
~ B(zi,t 2 ) ~ 



B(*2)<Tl(t 2 ) 



B(z 2 ,t 2 ) 



B(z p ,t 2 ) 



o-i(t 2 ). 



A natural question arises, when one can reconstruct S(X,t 2 ) € I, when one knows its pole data, i.e., if 
one knows all pole triples for all z's. 



8.3 Realization theorem (of Mittag-Leffler type) 

Let S(X,t 2 ) G I. Suppose that in order to solve a Mittag Leffler type problem, we are given a pole 
triple for each t 2 , (X(t 2 ), T, Y(t 2 )), where we have denoted the matrix T independent of t 2 . As we have 
seen, the pole pairs can be chosen so that the differential equations of lemma 18.21 with spectral matrix 
parameter T are satisfied. Then the following theorem answers the question of reconstructing S(X,t 2 ) 
from its pole triple. 

Theorem 8.4 Suppose that we are given a pole triple (C(i 2 ), A, B(t 2 )o~i) for each t 2 with constant 
A and C(t 2 ), B{t 2 )a\ solutions of the output VS. 8\) and the inverse input \5.3\) ODEs with the matrix 
spectral parameter A. Suppose also that an analytic function D(t 2 ) (the value at infinity) is given. Then 
the matrix function 

S(X, t 2 ) = D(t 2 ) + C(t 2 )(XI - A^B^ai (8.5) 

maps solutions of {3. 7y with spectral parameter X to solutions of US. 8\) with the same spectral parameter 
iff the linkage conditions i2.10.Linify are satisfied. 

Proof: The proof is essentially the same as for theorem 18.11 In order to obtain that the function 
S(X,t 2 ) built in (|8.5| maps solutions of (13. 7p to solutions of ()3.8|) with the same spectral parameter, it 
is necessary to prove that S(X,t 2 ) satisfies the following differential equation: 

J-S(A, t 2 ) = *^(* 2 *X + 7 *)S(A, h) - S(X, t 2 )a^(a 2 X + 7 ), (8.6) 
ot 2 

i.e., is of the form defined by (|3.9|) . Differentiating (|8.5p we obtain 

£- 2 S(X,t 2 ) = £- 2 D(t 2 ) + <r£fa.A + 7 ,)C(t 2 )(A/ - A^Bit,)^- 

-C{t 2 ){XI - A)- l B{t 2 ){a 2 A + 7 ). 

This can be rewritten as 

£- 2 S(X,t 2 ) = £- 2 D(t 2 ) -a^(t 2 )a 2 ,C(t 2 )B(t 2 )ax + C(t 2 )B(t 2 )a 2 (t 2 ) + a u 1 (t 2 )^(t 2 )D(t 2 )+ 
+X(a^(t 2 )a 2 *D(t 2 ) - Da^a 2 ) + a^(a 2 *X + ^){D{t 2 ) + C(t 2 ){XI - A) -1 B(t 2 )ai(t 2 ))- 

-(D(t 2 ) + C(t 2 )(XI - A)- 1 B(t 2 ))( ( r 2 A + 7) = 
= §- 2 D(t 2 ) - o-^{t 2 )o- 2 *C{t 2 )B(t 2 )o-i +C(t 2 )B(t 2 )a 2 (t 2 ) + o-^{t 2 ) 1 *{t 2 )D{t 2 )+ 
+X(a^(t 2 )or 2 *D(t 2 ) - Da^a 2 ) + X(a^(t 2 )a 2 ,S(X,t 2 ) - S(X,t 2 )a^(a 2 X + 7 ). 
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Demanding further that the differential equation (|8.6p holds for S(\,t 2 ), we obtain, considering first big 
A and then arbitrary one 

1- -a^a 2 *D(t 2 ) + D{t 2 )a^a 2 = and 

2. ^D(t 2 ) - a^D(t 2 ) + D(t 2 )a^ + a^a 2 *C{t 2 )B{t 2 )a x - C(t 2 )B(t 2 )a 2 = 0. 
Denoting D(t 2 ) = ai^D^^a^ 1 , we obtain that these conditions are 
f . -a 2 *D(t 2 ) + D{t 2 )a 2 = and 

2. ^D(t 2 ) - <J^^D(t 2 ) + trf^-D^fr + a^<J 2 *C{t 2 )B(t 2 )ax - C{t 2 )B{t 2 )a 2 - 0, 
which are exactly the linkage conditions (|2.10.Link[) . Thus the theorem holds. □ 
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